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Executive  Summary 


—  As  a  consequence  of  the  abbreviated  grant  period  (1 1  months),  we  did  not  make  as  much 
progress  as  we  anticipated  in  our  original  3-year  proposal. 

-  #  We  finished  and  published  a  study  on  c -fos  expression  in  the  SCN  of  tau  mutant 
hamsters.  We  demonstrated  a  correlation  between  phase-shifting  and  Fos  induction  thresholds  under 
conditions  where  both  responses  are  dramatically  altered  by  the  previous  light  history,  suggesting 
a  causal  association  between  changes  in  behavioral  phase  shifting  and  c -fos  in  the  SCN. 


•  We  finished  and  published  a  lengthy  study  of  multiunit  electrical  activity  (MU  A)  rhythms 
recorded  from  the  SCN  and  other  brain  areas  of  awake,  behaving  hamsters.  This  study  was 
something  of  a  technical  tour-de-force.  There  are  only  one  or  two  others  like  it  in  the  circadian 
literature](none  as  thorough).  The  conclusions  were  as  follows: 

o  Circadian  MUA  rhythms  in  the  SCN  had  periods  that  reflected  the  genotype  of  the 


subject. 

o  Circadian  MUA  rhythms  in  most  other  brain  areas  were  180°  out  of  phase  with 


those  in  the  SCN. 

o  Circadian  MUA  rhythms  in  the  bed  nucleus  were  in  phase  with  those  in  the  SCN. 
o  Ultradian  rhythms  with  2  periods  were  measured  (80  minutes  &14  minutes).  In 
contrast  to  the  circadian  rhythms,  these  were  not  affected  by  the  tau  mutation. 


—  •  We  modeled  the  pacemaker  system  of  the  tau  mutant  hamster  using  computer  simulation. 
These  simulations  led  to  two  primary  conjectures:  (1)  the  total  amplitude  of  the  pacemaker  system 
in  tau  mutant  hamsters  is  less  than  in  wild  type  animals,  and  (2)  the  coupling  between  the  unit  E  and 
M  oscillators  is  weakened  during  continuous  exposure  of  hamsters  to  DD.  * 
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Circadian  behavior  and  c-fos  expression 

In  hamsters  homozygous  for  the  circadian  clock  mutation  tau,  the  photic  history 
dramatically  affects  the  magnitude  of  light-induced  circadian  phase  shifts.  The  maximum 
amplitude  of  phase  shifts  produced  by  one-hour  light  pulses  presented  at  circadian  time  14  was 
less  than  2  hours  in  animals  that  had  been  in  constant  darkness  (DD)  for  2  days,  while  animals 
that  had  been  kept  in  DD  for  49  days  could  be  shifted  by  more  than  8  hours. 

In  this  study,  we  compared  the  effect  of  previous  light  history  on  the  amplitude  of 
circadian  phase  shifts  and  on  c-fos  expression  in  the  suprachiasmatic  nucleus  (SCN)  of  tau  mutant 
hamsters.  Although  the  maximum  amplitude  of  behavioral  phase  shifts  was  drastically  different 
between  animals  that  had  been  held  for  either  2  or  49  days  in  DD,  maximal  fos  induction  was  not 
significantly  different  in  these  two  groups.  However,  photic  thresholds  for  light-induced 
behavioral  phase  shifts,  c-fos  mRNA,  and  Fos  immunoreactivity  were  closely  correlated  within 
both  groups,  and  these  thresholds  were  lower  (more  sensitive  to  light)  after  49  than  after  2  days 
in  DD.  The  correlation  between  phase  shifting  and  Fos  induction  thresholds,  under  conditions 
where  both  responses  are  dramatically  altered  by  the  previous  light  history,  demonstrates  an 
association  between  changes  in  circadian  behavioral  phase  shifting  responses  of  tau  mutant 
hamsters  and  plasticity  of  light-induced  c-fos  expression  in  suprachiasmatic  nucleus.  However, 
because  the  maximum  amplitudes  of  Fos  induction  and  phase  shifting  were  not  correlated  in 
animals  that  had  been  in  DD  for  2  days,  we  speculate  that  the  level  of  c-fos  expression  does  not 
directly  determine  phase  shift  amplitude. 

In  vivo  recording  of  rhythmic  multiunit  activity 

We  recorded  multiple  unit  neural  activity  [multiunit  activity  (MU A)]  from  inside  and  outside 
of  the  suprachiasmatic  nucleus  (SCN)  in  freely  moving  male  golden  hamsters  housed  in  running- 
wheel  cages  under  both  light. /dark  cycles  and  constant  darkness.  The  circadian  period  of  MU  A  in 
the  SCN  matched  the  period  of  locomotor  activity:  it  was  ~24  hr  in  wild-type  and  20  hr  in 
homozygous  tau  mutant  hamsters.  The  peak  of  MU  A  in  the  SCN  always  occurred  in  the  middle 
of  the  day  or,  in  constant  darkness,  the  subjective  day.  There  were  circadian  rhythms  of  MUA 
outside  of  the  SCN  in  the  ventrolateral  thalamic  nucleus,  the  caudate  putamen,  the  accumbens 
nucleus,  the  medial  septum,  the  lateral  septum,  the  ventromedial  hypothalamic  nucleus,  the  medial 
preoptic  regions,  and  the  stria  medullaris.  The  rhythms  were  out  of  phase  with  the  electrical  rhythm 
in  the  SCN  but  in  phase  with  the  rhythm  of  locomotor  activity,  peaking  during  the  night  or 
subjective  night.  In  addition  to  circadian  rhythms,  there  were  significant  ultradian  rhythms  present; 
one,  with  a  period  of  ~80  min,  was  in  antiphase  between  the  SCN  and  other  brain  areas,  and  another, 
with  a  period  of  ~  1 4  min,  was  in  phase  between  the  SCN  and  other  brain  areas.  The  periods  of  these 
ultradian  rhythms  were  not  significantly  different  in  wild-type  and  tau  mutant  hamsters.  Of 
particular  interest  was  the  unique  phase  relationship  between  the  MUA  in  the  bed  nucleus  of  the  stria 
terminalis  (BNST)  and  the  SCN:  in  these  two  areas,  both  circadian  and  ultradian  components  were 
always  in  phase.  This  suggests  that  the  BNST  is  strongly  coupled  to  the  SCN  and  may  be  one  of  its 
major  output  pathways.  In  addition  to  circadian  and  ultradian  rhythms  of  MUA,  neural  activity  both 
within  and  outside  the  SCN  was  acutely  affected  by  locomotor  activity.  Whenever  a  hamster  ran 
on  its  wheel,  MUA  in  the  SCN  and  the  BNST  was  suppressed,  and  MUA  in  other  areas  was 
enhanced. 
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Modeling  the  pacemaker  of  the  tau  mutant  hamster 

Circadian  pacemakers  in  many  animals  are  compound.  In  rodents,  a  two-oscillator  model 
of  the  pacemaker  comprising  an  evening  ( E)  and  a  morning  (M)  oscillator  has  been  proposed  based 
on  the  phenomenon  of  “splitting”  and  bimodal  activity  peaks.  For  our  computer  simulations,  we 
viewed  the  hamster  pacemaker  as  a  system  of  mutually  coupled  E  and  M  oscillators.  Tau  mutant 
animals  exhibit  normal  type-1  phase-response  curves  (PRCs)  when  released  into  DD,  but  make  a 
transition  to  a  type-0  PRC  when  held  for  many  weeks  in  DD.  The  two-oscillator  model  describes 
particularly  well  some  behavioral  experiments  on  these  hamsters.  We  tried  to  determine  the 
relationships  between  oscillator  amplitude,  period,  PRC,  and  activity  duration  through  computer 
simulations.  Two  complementary  approaches  proved  useful  for  analyzing  weakly  coupled  oscillator 
systems.  We  adopted  a  “distinct  oscillators”  view  when  considering  the  component  E  and  M 
oscillators  and  a  “system”  view  when  considering  the  system  as  a  whole.  For  strongly  coupled 
systems,  only  the  system  view  is  appropriate.  The  simulations  lead  us  to  two  primary  conjectures: 
(1)  the  total  amplitude  of  the  pacemaker  system  in  tau  mutant  hamsters  is  less  than  in  the  wild-type 
animals,  and  (2)  the  coupling  between  the  unit  E  and  M  oscillators  is  weakened  during  continuous 
exposure  of  hamsters  to  DD.  As  coupling  strength  decreases,  activity  duration  (a)  increases  due  to 
a  greater  phase  difference  between  E  and  M.  At  the  same  time,  the  total  amplitude  of  the  system 
decreases,  causing  an  increase  in  observable  PRC  amplitudes.  Reduced  coupling  also  increases  the 
relative  autonomy  of  the  unit  oscillators.  The  relatively  autonomous  phase  shifts  of  E  and  M 
oscillators  can  account  for  both  immediate  compression  and  expansion  of  activity  bands  in  tau 
mutant  and  wild-type  hamsters  subjected  to  light  pulses. 


4 


Circadian  Behavior  and  Plasticity  of 
Light-Induced  c-fos  Expression  in 
SCN  of  tau  Mutant  Hamsters 

WoShtoomWUonM^ornh^ 

^T^terfvZosical ' 

fNSF  Center  for  IL  6°2°8' 

Physiology,  Northwes  wiLa!  Institute,  . 


tN  SF  Center  for ^^'^ersity,  E^nston,  IL  602U»,  t 

^■XZTh^ 

NSF^rfc^^ 


mr  Department 

"Current  address:  Sleep  Research  ‘ '  ce  Di„jLn,  National  Institute 

C“™  Current  addnes^ ^Zt^ha,  lharaht,  ,apan 

Abslmt  ta  hamsters  homozygous  £  £ ;  — ‘  ^ %£?££ 

photic  history  dramatically  af““  ,  ^ase  shifts  produced  by  1*¥< 

phase  shifts.  The  maximum  ampktu  P  ^  ^  had  been  m  DD  for  2 
pulses  presented  at  CT 14  was  less  than  2  h  ^  ^  days  could  be  shifted  by 

days,  whereas  animals  tkathadbeen  ep  ared  ^  e£fect  pf  previous  kght 

more  than  8  h.  In  this  study,  the  s^ts  and  on  c-fos  expression  m  the 

history  on  the  amplitude  of  maximum  amplitude  of  behaviora 

SCN  of  tau  mutant  hamsters.  Although  ^  that  had  been  held  for 

phase  shifts  was  drastically  difeent  was  not  significantly  different 

either  2  or  49  days  in  DD,  maxun^tic  thresholds  for  light-induced  behavioral 
in  these  two  groups.  Howe^h  ^f^oreactivity  were  closely  corrected 
phase  shifts,  c-fos  mRNA,  and  F  lower  (more  sensitive  to  hght) 

within  both  groups,  and  ^^  ^^pn  between  phase  shifting  and  Fos 
after  49  than  after  2  days  in  DD.  The  co  responses  are  dramatically 

S^n.hmsbolds.unde^onsw^^  associatio„  = 
altered  by  the  previous  hght ,  ^o^'  onses  of  tau  mutant  hamsters 

changes  l  circadian  behavioral  Ph^*£ss£  n  ?  SCN.  However,  because  die 
a  nlasticitv  of  light-induced  c  fos  p  chifting  were  not  correlated 

o,  POS ^ -  >eve,  of  <*. 

•  c-fbs,  lea  mulaut,  phase  shift,  light,  circadian,  SCN 

1  To  whom  all  cotiespondence  should  be  addressed. 

JOURNAL  OF  WOLOGICAL  RHYTHMS,  W>1- 13  No.  4,  Au^..«  »«• 


© 


1998  Sage  Publications,  Inc. 


306  JOURNAL  OF  BIOLOGICAL  RHYTHMS  /  August  1998 


INTRODUCTION 

A  variety  of  evidence  indicates  that  circadian  pace¬ 
makers  are  located  within  the  SCN  of  the  hypothala¬ 
mus  in  mammals  (Moore,  1995).  Entrainment  of  cir¬ 
cadian  rhythms  by  light  in  mammals  is  mediated  by 
retinal  photoreceptors  that  project  to  the  SCN  via  a 
direct  retinohypothalamic  tract.  Recent  studies  have 
shown  that  the  proto-oncogenes  c -/os  and  jun-B  are 
induced  in  the  rodent  SCN  in  response  to  light  (Aronin 
et  al.,  1990;  Earnest  et  al.,  1990;  Komhauser  et  al.,  1990, 
1992;  Rea,  1989;  Rusak  et  al.,  1990, 1992;  Takeuchi  et  al., 
1993).  Fos  and  Jun  family  proteins  appear  to  function 
in  cellular  signal  transduction  by  coupling  transient 
stimuli  to  gene  regulation  in  the  nucleus  (Muller  et  al., 
1984;  Greenberg  and  Ziff,  1984;  Lau  and  Nathans, 

1987) .  The  Fos  protein  dimerizes  with  Jim  family  pro¬ 
teins  to  form  a  complex,  AP-1  (Lee  et  al.,  1987; 
Bohmann  et  al.,  1987;  Rauscher  et  al.,  1988;  Curran  and 
Franza,  1988;  Sassone-Corsi  et  al.,  1988),  which  regu¬ 
lates  the  transcription  of  target  genes  (Curran  et  al., 

1988) .  Castel  et  al.  (1997)  have  recently  estimated  that 
light  pulses  which  produced  saturating  phase  shifts  of 
the  circadian  locomotor  rhythm  of  Sabra  mice  induce 
c-fos  expression  in  a  heterogeneous  population  of  SCN 
cells  comprising  about  one  fifth  of  the  total  number  of 
cells  in  the  SCN. 

The  characteristics  of  the  photic  induction  of  c-/os 
and  jun- B  in  the  SCN  are  correlated  in  several  respects 
with  behavioral  phase-shifting  responses  (Komhauser 
et  al,  1996).  First,  the  induction  of  c-fos  and  jun- B  is 
limited  to  the  retinorecipient  area  of  the  SCN.  Second, 
the  light  sensitivities  of  the  behavioral  phase-shifting 
response  and  of  the  induction  of  c-fos  mRNA  are  quan¬ 
titatively  correlated.  Third,  c-fos  and/wn-B  are  induced 
by  light  only  during  the  subjective  night,  a  time  when 
light  pulses  cause  phase  shifts  in  circadian  rhythms. 
Fourth,  injection  of  antisense  oligonucleotides  into  the 
third  ventricle  (which  block  the  translation  of  c-fos  and 
jun- B  in  the  SCN)  before  a  light  pulse  prevents  light- 
induced  phase  delays  in  rats.  These  results  suggest 
that  the  light-induced  expression  of  these  two  genes 
may  be  causally  involved  in  photic  entrainment  of 
circadian  rhythms. 

The  circadian  period  mutation,  tau,  affects  both  the 
free-running  period  and  the  shape  of  the  phase  re¬ 
sponse  curve  in  the  golden  hamster  (Ralph  and 
Menaker,  1988;  Shimomura  and  Menaker,  1994).  tau 
mutant  hamsters  regularly  exhibit  high-amplitude 
phase  resetting  (type  0  resetting)  (Winfree,  1980)  under 
certain  circumstances,  but  it  is  rare  to  find  type  0  reset¬ 


ting  in  wild  type  hamsters.  The  magnitude  of  light- 
induced  phase  shifts  in  homozygous  tau  mutant  ham¬ 
sters  is  dramatically  affected  by  the  animal's  light 
history  prior  to  receiving  a  light  pulse.  For  example, 
after  short  durations  (2  or  7  days)  in  DD,  the  phase 
response  curve  is  type  1  (maximum  amplitude  about 
2-3  h),  phase  resetting  does  not  resemble  type  0,  even 
in  response  to  a  bright  light  pulse,  whereas  after  long 
durations  in  DD  (e.g.,  49  days),  it  becomes  type  0 
(maximum  amplitude  about  8-10  h).  Grosse  et  al. 
(1995)  found  that  although  the  magnitude  of  phase 
advances  produced  by  saturating  light  pulses  was 
much  larger  in  homozygous  tau  mutant  hamsters  than 
wild  types,  the  number  of  Fos  immunoreactive  (IR) 
cells  in  the  SCN  was  similar  between  the  two  geno¬ 
types.  Their  study  suggests  that  Fos  may  not  deter¬ 
mine  the  amplitude  of  light-induced  phase  shifts 
directly 

We  report  here  the  results  of  experiments  designed 
to  determine  whether  change  in  circadian  behavioral 
responses  between  2  and  49  days  in  DD  in  tau  mutant 
hamsters  is  accompanied  by  differences  in  photic  sen¬ 
sitivity  and,  if  so,  whether  the  differences  in  both 
sensitivity  and  maximum  response  are  correlated  with 
c-fos  expression  in  the  SCN. 

METHODS 
Animals  and  Housing 

Homozygous  tau  mutant  hamsters  were  bred  from 
homozygous  parents  at  the  University  of  Virginia.  All 
animals  in  these  experiments  were  raised  on  a  20-h  LD 
cycle  (L:D  =  11.7:8.3,  which  has  the  same  L:D  ratio  as 
L:D  =  14:10)  since  birth.  After  weaning,  they  were 
group  housed  (1-3  animals  per  cage).  At  5-12  weeks  of 
age,  they  were  transferred  to  individual  cages  (20  x  20 
x  47  cm)  equipped  with  running  wheels  (8  cm  in 
width,  17  cm  in  diameter).  Running- wheel  cages  hous¬ 
ing  single  animals  were  kept  in  ventilated,  light-tight 
boxes  (six  cages  per  box),  initially  under  the  same  LD 
cycles  in  which  the  animals  were  raised.  Illumination 
(provided  by  fluorescent  lights)  was  approximately 
10-40  jiW  cm“2  at  15  cm  above  the  cage  bottom.  Food 
and  water  were  available  ad  libitum.  Ambient  tem¬ 
perature  was  approximately  21°C.  Cages  were 
changed  at  least  every  28  days,  and  water  bottles  every 
14  days,  without  using  visible  light  (infrared  viewer, 
FJW  Optical  Systems,  Palatine,  IL,  USA)  for  hamsters 
held  in  DD.  Each  hamster  was  used  only  once.  All 
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Figure  1.  Locomotor  activity  records  from  homozygous  tau  mu¬ 
tant  hamsters  kept  for  7  days  (A  &  C)  or  49  days  (B  &  D)  in  DD. 
The  animals  whose  records  are  shown  in  A  and  B  received  a  light 
pulse  of  about  200  cm"2  at  CT  14;  those  whose  records  are 
shown  in  C  and  D  received  the  same  irradiance  light  pulse  at  CT 
18.  The  times  of  the  light  pulses  are  indicated  by  asterisks;  the  data 
are  plotted  on  a  20-h  time  base,  and  only  the  7  days  before  and  the 
14  days  after  the  light  pulse  are  shown. 

experiments  were  performed  in  accordance  with  the 
Institute  of  Laboratory  Animal  Resources  guidelines 
set  forth  in  the  1985  publication  Guide  for  the  Care  and 
Use  of  Laboratory  Animals . 

Recording  of  Locomotor  Activity 

Wheel-running  activity  was  recorded  with  micro 
switches  activated  by  the  wheel  (using  Data  Quest  III, 
Data  Sciences  International,  St.  Paul,  MN,  USA);  activ¬ 
ity  was  counted  as  the  number  of  wheel  revolutions 
per  6  min. 

Photic  Stimulation  and  Circadian  Phase 

Throughout  this  manuscript,  the  word  "days"  is 
used  to  refer  to  circadian  cycles  (whether  approxi¬ 
mately  20  h  or  24  h).  After  at  least  10  days  of  steady 
state  entrainment,  hamsters  were  released  into  DD. 
The  number  of  days  in  DD  was  varied  systematically, 
and  a  light  pulse  of  a  specific  irradiance  was  then 
presented  at  an  appropriate  phase  of  the  circadian 
cycle.  Light  bulbs  were  wrapped  with  black  plastic 
when  reduced  irradiance  was  required.  When  it  was 


Figure  2.  Relationship  between  magnitudes  of  phase  shifts  (open 
bars)  and  numbers  of  Fos  immunoreactive  (IR)  cells  in  the  SCN 
(shaded  bars).  Phase  shifts  were  induced  by  a  60-min  bright  (2.0 
x  102  pW  cm”2)  light  pulse  in  tau  mutant  hamsters  after  2  days  (A) 
or  49  days  (B)  in  DD.  Fos  immunoreactivity  was  examined  at  the 
end  of  the  60-min  light  pulse.  Each  data  point  represents  the  mean 
± SEM  of  5-6  animals  for  phase  shifts  and  3-4  animals  for  Fos  IR. 
Representative  Fos  immunostaining  in  the  SCN  is  shown  for 
animals  given  a  light  pulse  after  2  days  (C)  or  49  days  (D)  in  DD. 


necessary  to  give  animals  light  stimuli  of  different 
irradiances;  animals  were  put  in  small  plastic  cups  (10 
cm  diameter,  20  cm  height),  which  were  moved  to  a 
light-tight  box  illuminated  by  cool  white  fluorescent 
light  for  60  min;  otherwise  they  were  left  in  their  own 
cages  and  the  cages  were  moved  into  the  same  type  of 
light-tight  box.  Light  irradiance  was  measured  15  cm 
above  the  bottom  of  the  cup  or  cage.  By  convention, 
each  day  of  a  free-running  circadian  rhythm  is  divided 
into  24  h  of  CT.  Therefore,  the  duration  of  1  circadian 
hour  is  obtained  by  dividing  the  period  of  the  cir¬ 
cadian  rhythm  of  an  individual  animal  by  24  (e.g.,  if 
the  free-running  period  is  20  clock  hours,  1  circadian 
hour  is  50  clock  minutes).  In  nocturnal  animals,  activ¬ 
ity  onset  is  defined  as  CT  12  (thus,  CT  18  occurs  6 
circadian  hours  after  activity  onset).  The  magnitude 
and  direction  (delay  or  advance)  of  phase  shifts  de¬ 
pend  upon  the  circadian  phase  of  stimulation.  Maxi¬ 
mum  phase  delays  are  produced  by  light  stimuli  at  CT 
13-15  and  maximum  phase  advances  are  produced 
at  CT  18-20  in  homozygous  tau  mutant  hamsters 
(Shimomura  and  Menaker,  1994).  In  this  study,  light 
stimuli  beginning  at  CT  14  were  used  to  produce 
phase  delays  and  light  stimuli  beginning  at  CT  18  to 
produce  phase  advances.  All  light  stimuli  were  1  h  in 
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Figure  3.  Tune  course  of  Fos  inununoreactive  (IR)  in  the  SCN  of 
tau  mutant  hamsters  following  a  light  pulse.  Animals  received  a 
60-min  light  pulse  (of  2.0  x  102  pW  cm'2  irradiance)  at  CT 14.  Open 
circles  represent  Fos  IR  in  animals  that  had  been  maintained  for  2 
days  in  DD,  and  closed  circles  represent  animals  previously  kept 
for  49  days  in  DD.  Each  data  point  represents  the  mean  ±SEM  of 
the  number  of  Fos  IR  cells  in  the  SCN  of  three  animals.  The  x-axis 
indicates  time  after  the  start  of  the  60-min  light  pulse. 


duration.  After  light  stimulation,  hamsters  were  re¬ 
turned  (using  an  infrared  viewer)  to  DD  for  at  least  14 
days  (for  measurement  of  phase  shifts)  or  killed  in 
darkness  (for  immunocytochemistry  or  in  situ  hy¬ 
bridization).  During  this  period,  cages  were  not 
changed,  so  as  to  prevent  phase  shifts  (which  can  be 
produced  by  cage  changes)  (Reebs  and  Mrosovsky, 
1989). 

Data  Analysis 

The  magnitude  of  phase  shifts  was  determined  by 
measuring  the  phase  difference  (using  activity  onset 
as  a  phase  reference  point)  between  eye-fitted  regres¬ 
sion  lines  through  3-7  consecutive  activity  onset  times 
immediately  prior  to  the  light  pulse,  and  at  least  10 
consecutive  activity  onset  times  after  the  light  pulse 
(excluding  the  four  cycles  immediately  after  the 
pulse).  The  magnitude  of  all  phase  shifts  is  expressed 
in  circadian  hours.  Differences  between  two  groups 
were  evaluated  by  Student’s  f-test.  Significant  differ¬ 
ence  was  defined  as  p  <  .05. 

Immunocytochemistry 

Animals  were  anesthetized  with  a  lethal  dose  of 
chloroform  by  inhalation  and  then  perfused  intra- 


cardially,  first  with  20-30  ml  physiological  saline 
(0.9%)  containing  150 IU  heparin/ 10  ml,  followed  by 
250  ml  2%  formaldehyde  in  0.01  M  phosphate  buffer 
saline  (PBS;  pH  7.4)  with  picric  acid  (15%).  After  per¬ 
fusion,  brains  were  removed  and  post  fixed  at  4°C. 
Serial  coronal  sections  were  cut  at  30  pm  or  50  pm 
using  a  cryostat  and  collected  m  PBS.  Free-floating 
sections  were  transferred  to  a  solution  of  normal  goat 
serum  (1  x  30'1  in  PBS-T)  for  30-60  min  then  directly 
into  primary  antiserum.  The  Fos  antiserum  used  in 
this  study  was  an  anti-Fos  (4-17)  rabbit  polyclonal 
antiserum  purchased  from  Oncogene  Science 
(Manhasset,  NY,  USA).  The  antibodies  were  used  at 
a  dilution  of  1:2000  for  48-72  h  at  4°C.  The  sites  of 
antigen-antibody  binding  were  visualized  with  an 
avidin-biotin-peroxidase  procedure  (Elite  ABC  kit, 
Vector  Labs,  Burlingame,  CA,  USA).  Sections  were 
washed  in  Tris  buffer  (pH  7.4)  before  incubation  for  3 
min  in  0.025%  diaminobenzidine  (DAB)  containing 
0.003%  peroxide.  Sections  mounted  on  gelatin-coated 
slides  were  dehydrated  through  graded  alcohols  into 
xylene  and  left  overnight  before  being  rehydrated  to 
distilled  water  and  immersed  in  0.2%  osmium  tetrox- 
ide  solution  for  2  min  to  intensify  the  DAB  reaction. 
Labeled  cell  nuclei  in  the  ventrolateral  subdivisions  of 
SCN  in  six  consecutive  sections  were  counted  without 
knowledge  of  the  experimental  conditions  to  which 
the  animal-had  been  exposed. 

In  Situ  Hybridization 

In  situ  hybridization  was  performed  essentially  as 
previously  described  (Komhauser  et  al.,  1990).  Briefly, 
30  min  after  the  beginning  of  the  light  pulse,  animals 
were  killed  by  decapitation  and  their  eyes  were  re¬ 
moved  in  darkness.  Brains  were  quickly  removed  in 
dim  red  light,  frozen  on  dry  ice,  and  stored  at  -80°C 
until  sectioning.  Twenty  pm-thick  coronal  sections 
were  cut  and  mounted  on  gelatin/poly-L-lysine 
coated  slides.  Sections  were  air  dried,  fixed  in  5% 
paraformaldehyde,  and  hybridized  for  18  h  at  47°C  in 
a  solution  of  50%  formamide,  300  mM  NaCl,  10  mM 
Tris-HCl  (pH  8.0),  1  mM  EDTA,  IX  Denhardt's  solu¬ 
tion,  10%  dextran  sulfate,  10  mM  DTT,  0.5  pg  ml'1 
tRNA,  0.5  pg  ml'1  poly(A),  and  2xl07  cpm  ml-1  of 
35S-labeled  c -fos  riboprobe  (transcribed  from  a  full- 
length  mouse  c-fos  cDNA  clone  in  a  pGEM-1  vector, 
linearized  with  Bgl  II  to  produce  an  anti-sense  tem¬ 
plate  of  approximately  1.8  kb  in  length,  or  with  Nco  I 
for  a  sense-strand  template  of  approximately  1.1  kb). 
Following  treatment  with  20  pg  ml'1  RNase  A  in  2X 
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Figure  4.  Tune  course  of  c -/os  mRNA  in  the  SCN  of  tau  mutant 
hamsters  after  a  light  pulse.  Animals  received  a  60-min  light  pulse 
(200  LiW  cm~2  irradiance)  at  CT 14.  Open  circles  represent  relative 
mRNA  levels  in  animals  previously  maintained  for  2  days  in  DD, 
and  closed  circles  represent  the  levels  in  animals  previously  main¬ 
tained  for  49  days  in  DD.  Each  point  represents  the  mean  ±SEM 
of  c-fos  mRNA  levels  determined  by  in  situ  hybridization  (plotted 
relative  to  the  maximum  value).  The  x-axis  indicates  time  after  the 
start  of  the  60-min  light  pulse. 


Log  Irradiance  (jiw/cm  ) 


SSC  (Standard  Saline  Citrate)  at  37°C  for  30  min,  slides 
were  washed  successively  in  2X  SSC,  IX  SSC,  0.5X 
SSC,  and  0.1X  SSC  for  30  min  each  at  4 7°C  and  were 
then  dehydrated  through  an  ethanol  series.  Dried 
slides  were  first  exposed  to  X-ray  film  (Amersham  Hy¬ 
perfilm  b-max),  and  then  to  autoradiographic  emul¬ 
sion  (Kodak  NTB-2)  for  3  weeks.  After  developing, 
sections  were  stained  with  cresyl  violet. 

Quantification  of  specific  hybridization  was  per¬ 
formed  using  the  Image-1/  AT  program  from  Univer¬ 
sal  Imaging  Corporation,  as  described  previously 
(Komhauser  et  al.,  1990).  From  each  brain,  all  intact 
sections  containing  portions  of  the  SCN  were  included 
in  our  analysis.  Using  100X  dark-field  magnification, 
an  area  of  75  x  75  pm  that  contained  the  highest  silver 
grain  density  in  each  SCN  was  analyzed  using  an 
" area-brightness''  measurement.  To  normalize  each 
SCN  measurement  for  variations  in  background,  an 
area  of  equal  size,  225  pm  lateral  to  the  first  (outside 
the  SCN)  was  also  measured.  We  defined  the  differ¬ 
ence  between  the  measurement  within  the  SCN  and 
that  in  the  lateral  area  as  the  "signa 1"  in  each  section. 
Values  representing  the  mean  signals,  in  both  SCNs, 
of  all  sections  from  one  brain  were  then  normalized  by 
taking  their  ratio  to  those  obtained  from  hamsters  that 
received  no  light. 


Figure  5.  The  magnitude  of  phase  shift  (A)  and  Fos  induction  (B) 
in  the  tau  mutant  SCN  as  a  function  of  light  pulse  irradiance. 
Animals  that  had  been  previously  held  for  either  2  or  49  days  in 
DD  received  60-min  light  pulses  of  varying  irradiance  at  CT  14. 
Animals  used  for  measuring  phase  shifts  were  kept  in  running 
wheels  for  14  days  in  DD  after  the  light  pulse,  those  for  immuno- 
cytochemistry  were  killed  at  the  end  of  the  60-min  light  pulse,  and 
those  for  in  situ  hybridization  were  killed  30  min  after  light  onset. 
In  panel  A,  open  drdes  represent  the  magnitude  of  phase  shifts 
after  2  days  in  DD  and  closed  circles  represent  the  magnitude  of 
phase  shifts  after  49  days  in  DD.  Each  data  point  represents  the 
mean  ±SEM  of  5-6  animals.  In  panel  B,  open  circles  represent  the 
number  of  Fos  immunoreactive  (IR)  cells  in  animals  previously 
kept  in  DD  for  2  days,  closed  circles  represent  the  number  of  Fos 
IR  cells  in  animals  previously  kept  in  DD  for  49  days,  open 
triangles  represent  relative  c-fos  mRNA  in  animals  previously 
kept  in  DD  for  2  days,  and  closed  triangles  represent  relative  c -fos 
mRNA  in  animals  previously  kept  in  DD  for  49  days.  Each  data 
point  represents  the  mean  ±SEM  of  3-4  animals. 

RESULTS 

Magnitude  of  Phase  Shift 

To  determine  whether  changes  in  the  magnitude  of 
behaviorally  measured  phase  shifts  are  paralleled  by 
changes  in  the  expression  of  c-/os  in  the  SCN  of  tau 
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Figure  6.  Fos  immun  ore  active  (IR)  in  the  SCN  and  phase  shift 
responses  plotted  as  a  function  of  stimulus  irradiance  in  animals 
that  had  been  kept  in  DD  for  2  days  (A)  or  49  days  (B).  60-min  light 
pulses  of  varying  irradiance  were  given  at  CT  14,  and  Fos  IR  was 
determined  at  the  end  of  the  60-min  light  period.  Closed  circles 
represent  the  mean  ±SEM  of  numbers  of  Fos  IR  cells  from  3-4 
animals  in  each  group,  and  open  circles  represent  the  mean  (±S£Af) 
amplitude  of  light-induced  phase  shifts  from  5-6  animals  in  each 
group. 

mutant  hamsters,  c-fos  induction  was  examined  in 
animals  previously  held  for  2  days  or  49  days  in  DD. 
Fig.  1  shows  typical  light-induced  phase  shifts  in  tau 
mutant  hamsters.  After  2  days  in  DD  (Fig.  1  A,C),  the 
magnitude  of  phase  shifts  was  much  smaller  than  after 
49  days  in  DD  (Fig.  1  B,D),  at  both  CT  14  and  CT  18. 
Since  this  difference  was  larger  at  CT  14  (Fig.  1  A  vs. 
B),  we  used  light  pulses  at  CT  14  to  investigate  possible 
difference  in  c -fos  expression.  While  the  magnitude  of 
the  fnaximum  phase  shifts  induced  by  the  brightest 
light  pulses  (2.0  x  102  pW  cm"2)  was  dramatically  larger 
after  49  days  in  DD  than  after  2  days  in  DD,  the  number 
of  Fos  IR  cells  in  the  SCN  detected  at  60  min  after  the 
onset  of  light  was  not  significantly  different  between  the 
two  groups  (Fig.  2).  We  next  investigated  the  time 


course  of  Fos  immunoreactivity  to  determine  if  there 
were  differences  in  rates  of  synthesis  or  degradation 
between  the  two  groups  (Fig.  3).  After  induction  by 
light,  the  decline  in  amount  of  Fos  IR  appeared  to  be 
slightly  more  gradual  after  49  days  than  after  2  days  in 
DD.  Next  we  investigated  the  c-/os  mRNA  time  course 
after  2  days  or  49  days  in  DD.  The  mRNA  time  courses 
were  essentially  similar  in  the  two  groups  (Fig.  4). 

Threshold  for  Phase  Shift 

To  determine  the  effect  of  the  animal's  light  history 
on  the  photic  thresholds  for  phase  shifting,  Fos  IR,  and 
c-fos  mRNA  induction,  we  gave  light  pulses  of  varying 
irradiance  at  CT  14  to  animals  previously  held  in  DD 
for  2  or  49  days  (Figs.  5  and  6).  The  thresholds  for  phase 
shifting  (Fig.  5A)  and  for  induction  of  Fos  IR  and  c-fos 
mRNA  (Fig.  5B)  were  both  lower  after  49  days  than 
after  2  days  in  DD.  As  in  the  first  set  of  experiments, 
while  the  maximum  phase  shift  amplitude  was  obvi¬ 
ously  greater  after  49  days  in  DD,  the  maximum  num¬ 
ber  of  Fos-positive  cells  was  not  significantly  different. 

The  stimulus-response  curves  for  phase  shifting 
and  for  Fos  induction  are  compared  directly  in  Fig.  6, 
which  are  replotted  from  data  in  Fig.  5.  After  2  days  in 
DD,  the  effectiveness  of  light  for  phase  shifting  and 
for  Fos  induction  was  closely  correlated  for  light 
pulses  of  the  lowest  three  irradiance  levels,  but  this 
correlation  disappeared  at  the  two  highest  irradiances 
because  the  magnitude  of  the  phase  shift  response 
saturated  at  1.2  x  10°  p.W  cm'2  while  the  number  of  Fos 
IR-positive  cells  continued  to  increase  (Fig.  6A).  After 
49  days  in  DD,  the  thresholds  for  both  phase-shifting 
and  Fos-induction  responses  were  lower.  In  addition, 
a  significant  change  in  the  relationship  between  the 
two  responses  was  observed;  as  in  the  2-day  DD  ani¬ 
mals,  the  thresholds  for  phase  shifting  and  Fos-IR 
responses  were  indistinguishable;  however,  in  the  49- 
day  DD  animals  (but  not  in  the  2-day  DD  animals), 
both  responses  saturated  at  the  same  irradiance  (3.0  x 
1CT1  |iW  cm'2)  and  remained  closely  correlated  at 
higher  irradiances  (cf.  Fig.  6  A,B). 

DISCUSSION 

In  tau  mutant  hamsters,  the  magnitude  of  light- 
induced  circadian  phase  shifts,  as  well  as  the  irradi¬ 
ance  threshold  for  phase  shifting,  are  dramatically 
affected  by  the  prior  light  history,  especially  by  the 
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Figure  7.  Diagrammatic  illustration  of  the  changes  in  the  relationship  between  Fos  induction  and  phase  shift  responses  in  hamsters  held 
for  2  or  49  days  in  DD  (see  text  for  discussion). 


time  spent  in  DD  (Shimomura  and.  Menaker,  1994).  In 
this  further  study  of  the  tau  mutant  hamster  light 
response,  we  found  that  the  maximum  level  of  Fos 
induction  in  the  SCN  in  response  to  a  bright  light  pulse 
was  not  affected  by  the  duration  of  prior  exposure  to 
DD.  However,  the  sensitivity  of  the  light-induced  ex¬ 
pression  of  c-fos  in  the  SCN  was  dependent  upon  the 
duration  of  prior  exposure  to  DD.  The  irradiance 
threshold  for  c-/os  induction  was  much  lower  after  49 
days  than  after  2  days  in  DD.  This  decrease  in  thresh¬ 
old  (i.e.,  increase  in  light  sensitivity)  for  the  induction 
of  c-fos  mRNA  and  Fos  IR  paralleled  the  change  in 
threshold  observed  for  behavioral  phase  shifting. 

These  changes  in  the  relationship  between  Fos  in¬ 
duction  and  phase  shift  responses  are  illustrated  dia- 
grammatically  in  Fig.  7.  After  2  days  in  DD,  the  thresh¬ 
olds  for  phase  shifts  and  for  Fos  induction  are 
correlated,  but  the  phase  shift  response  reaches  satu¬ 
ration  with  increasing  irradiance,  whereas  Fos  induc¬ 
tion  continues  to  increase  up  to  the  highest  irradiance 
that  can  be  produced  by  our  experimental  apparatus. 
After  49  days  in  DD,  the  thresholds  for  both  responses 
are  lower  and  the  responses  increase  in  parallel.  How¬ 


ever,  while  -the  maximum  Fos  induction  is  approxi¬ 
mately  the  same  as  after  2  days  in  DD,  the  amplitudes 
of  phase  shifts  are  dramatically  larger.  Thus,  with 
prolonged  time  in  DD,  Fos  induction  shows  an  in- 
crease  in  sensitivity  but  the  same  maximum  response, 
while  the  phase  shift  response  increases  both  in  sensi¬ 
tivity  and  in  maximum  response. 

It  is  not  clear  whether  such  differences  in  both 
phase  shift  and  c-fos  expression  will  be  found  to  be 
general  characteristics  when  comparing  type  1  (low 
amplitude)  and  type  0  (high  amplitude)  phase  reset¬ 
ting.  There  are  several  ways  other  than  prolonged 
exposure  to  DD  to  produce  type  0  resetting  in  tau 
mutant  hamsters.  For  example,  type  0  resetting  can  be 
produced  by  pre-entrainment  to  some  but  not  other 
T-cycles.  When  animals  are  entrained  to  LD  =  1:19, 
phase  shifts  of  up  to  10  h  can  be  induced  by  light  pulses 
at  CT  14,  whereas  animals  pre-entrained  to  LD  =  1:18 
and  given  identical  light  pulses  show  shifts  of  less  than 
1  h  (Shimomura  and  Menaker,  1994).  Comparison  of 
c-fos  induction  in  these  two  groups  could  be  used  to 
distinguish  between  direct  effects  of  prolonged  dark¬ 
ness  and  the  effects  of  prolonged  darkness  on  cir- 
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cadian  organization.  Because  type  0  resetting  is  rarely 
observed  in  wild  type  hamsters,  it  would  also  be  in¬ 
formative  to  determine  whether  any  circumstances 
exist  in  which  the  phenomena  observed  in  this  study 
occur  in  wild  type  hamsters  (see  Elliott,  1994). 

The  induction  by  light  of  the  immediate  early  genes 
c-fos  and  jun-B  in  the  SCN  is  closely  correlated  in 
several  respects  with  light-induced  phase  shifting  of 
behavioral  circadian  rhythms  (for  review,  see  Kom- 
hauser  et  al.,  1996).  Because  of  this  relationship,  it  has 
been  proposed  that  transcription  factors  such  as  AP-1 
(composed  of  Fos/JunB  dimers)  could  play  a  role  in 
transducing  the  light  signal  into  changes  in  the  expres¬ 
sion  of  genes  involved  in  the  circadian  timing  mecha¬ 
nism  (Komhauser  et  al.,  1992;  Takeuchi  et  al.,  1993). 
Such  hypotheses  should  give  rise  to  specific  predic¬ 
tions  about  the  relationship  between  Fos  expression 
and  phase  shifting,  testable  in  a  variety  of  circum¬ 
stances  where  the  circadian  system  shows  distinctive 
or  atypical  responses  to  light.  In  the  current  study, 
dramatic  plasticity  is  observed  in  two  features  of  the 
phase-shifting  response— the  magnitude  of  phase 
shifts  and  light  sensitivity.  Therefore,  correlated 
changes  in  the  magnitude  and  sensitivity  of  Fos  induc¬ 
tion  can  be  proposed  as  criteria  for  testing  the  simplest 
possible  type  of  causal  relationship:  the  hypothesis 
that  Fos  induction  (acting  in  concert  with  other  tran¬ 
scription  factors)  results  in  phase  shifts  and  also  deter¬ 
mines  the  magnitude  of  those  shifts. 

Our  results  indicate  that  the  dramatic  increase  in 
the  amplitude  of  phase  shifts  observed  after  pro¬ 
longed  exposure  to  DD  is  not  paralleled  by  a  similar 
increase  in  Fos  induction.  The  number  of  Fos-positive 
cells  does  not  directly  correlate  with  the  magnitude  of 
light-induced  phase  shifts.  This  result  is  consistent 
with  previous  work  that  compared  the  number  of  Fos 
IR  cells  between  wild  and  tau  mutant  hamsters 
(Grosse  et  al.,  1995). 

On  the  other  hand,  we  found  that  even  though  the 
thresholds  for  induction  of  behavioral  phase  shifts 
and  c-fos  (both  mRNA  and  protein)  expression  were 
lower  after  49  days  than  after  2  days  in  DD,  the  thresh¬ 
olds  for  both  responses  were  closely  correlated  under 
all  conditions.  Komhauser  et  al.  (1990)  have  reported 
that  .the  thresholds  for  inducing  behavioral  phase 
shifting  and  c -fos  mRNA  are  quantitatively  correlated 
in  wild  type  hamsters.  The  present  study  demon¬ 
strates  that  in  tau  mutant  hamsters,  the  irradiance 
threshold  for  induction  of  both  c-fos  mRNA  and  Fos  IR 
in  the  SCN  is  correlated  with  that  for  induction  of 


phase  shifts,  and  furthermore  that  this  correlation  is 
maintained  when  the  phase-shifting  threshold  is 
changed  by  prior  history.  The  thresholds  for  inducing 
behavioral  phase  shifts  and  Fos  by  light  also  change 
during  aging  in  golden  hamsters  (Zhang  et  al.,  1996). 
Young  animals  (~8  weeks  of  age)  are  more  sensitive  to 
light  than  old  animals  (18  months  old)  in  both  behav¬ 
ioral  phase  shifting  and  Fos  induction  by  light.  In  that 
study  as  in  the  present  one,  even  when  the  threshold 
for  behavioral  phase  shifts  is  altered  by  the  animal's 
history,  the  threshold  for  c-fos  induction  is  changed  in 
the  same  direction  and  by  the  same  amount. 

These  two  sets  of  results  are  consistent  with  the  idea 
that  Fos  expression  may  be  involved  in  the  signaling 
mechanism  by  which  light  produces  phase  shifts. 
Ho-  ’ver,  they  indicate  that  if  Fos  does  participate  in 
the  .  nase-shifting  mechanism,  it  may  act  either  as  a 
trigger  or  as  a  permissive  signal,  but  that  the  level  of 
Fos  produced  does  not  by  itself  determine  the  magni¬ 
tude  of  the  phase  shift.  Stated  differently,  although 
Fos/JunB  dimers  may  act  to  regulate  gene  expression, 
which  ultimately  leads  to  phase  shifts,  the  amplitude 
of  the  response  to  this  input  pathway  is  likely  to  be 
regulated  by  some  other  (downstream)  component. 

The  idea  that  a  signaling  component  other  than  Fos 
normally  limits  (or  perhaps  even  regulates)  the  ampli¬ 
tude  of  phase  shifts  is  consistent  with  analysis  of  the 
irradiance  response  curves  obtained  in  the  present 
study  from  animals  that  had  been  exposed  to  different 
light  regimes  (Fig.  7).  In  animals  that  had  been  in  DD 
for  only  2  days,  the  amount  of  Fos  induction  increases 
continuously  with  increasing  irradiance,  far  beyond 
the  irradiance  at  which  the  phase-shifting  response 
reaches  saturation.  Therefore,  the  amplitude  of  the 
phase  shift  must  be  limited  by  something  other  than 
Fos.  In  contrast,  in  animals  that  had  been  in  DD  for  49 
days,  Fos  induction  attains  similar  levels  to  those  ob¬ 
served  in  2-day  DD  animals,  but  phase  shift  magni¬ 
tude  now  rises  in  concert  with  Fos  to  much  higher 
levels.  One  interpretation  of  this  difference  is  that  after 
prolonged  periods  in  DD,  some  change  occurs  such 
that  the  hypothetical  other  signaling  component  is  no 
longer  limiting,  and  higher  Fos  levels  can  produce  the 
observed  larger  phase  shifts.  This  other  component 
could  be  an  integral  part  of  the  circadian  timing 
mechanism  itself  or  part  of  the  light  input  pathwa-  o 
the  clock.  Mutations  such  as  tau  that  affect  circac  i 
period  are  likely  to  produce  alterations  in  the  funci  i 
of  one  or  more  molecular  components  of  the  met  - 
nism  that  generates  circadian  oscillations.  This  muca- 
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tion  also  changes  the  amplitude  of  the  PRC  to  non- 
photic  stimuli  (Mrosovsky  et  al.,  1992).  Period  muta¬ 
tions  that  are  known  to  affect  molecular  components 
of  the  clock  mechanism  in  both  Drosophila  and 
Neurospora  have  been  shown  to  affect  PRC  amplitude 
(Konopka  et  al.,  1989;  Dharmananda,  1980).  The  fact 
that  large  changes  in  the  amplitude  of  behavioral 
phase  shifts,  dependent  upon  light  history,  are  ob¬ 
served  more  generally  in  tau  mutant  than  in  wild  type 
•hamsters  (but  see  Elliott,  1994)  is  consistent  with  the 
suggestion  that  these  changes  are  a  consequence  of 
pacemaker  characteristics  that  have  been  altered  by 
the  mutation. 

Recently,  putative  mammalian  orthologs  to  the 
Drosophila  period  gene  have  been  reported  (called 
mPerl  and  mPer2)  (Sun  et  al.,  1997;  Tei  et  al.,  1997, 
Albrecht  et  al.,  1997;  Shearman  et  al.,  1997).  mPerl  and 
mPerl  mRNA  show  circadian  oscillations  in  mouse 
SCN,  and  mPerl  was  rapidly  induced  by  a  brief  light 
pulse  in  the  subjective  night  (Albrecht  et  al.,  1997; 
Shigeyoshi  et  al.,  1997;  Shearman  et  al.,  1997).  The 
photic  threshold  for  behavioral  phase  shifting  and 
mPerl  expression  in  the  SCN  were  quantitatively  cor¬ 
related  (Shigeyoshi  et  al.,  1997),  and  the  time  course  of 
mPerl  expression  was  slightly  slower  than  that  of  c-fos 
expression  (Albrecht  et  al.,  1997).  Therefore,  it  will  be 
of  great  interest  to  compare  mPerl  induction  by  light 
in  the  SCN  of  tau  mutant  hamsters  that  have  been  held 
in  DD  for  2  or  49  days  as  a  test  of  the  possibility  that 
■mPerl  expression  is  being  influenced  by  immediate 
early  genes. 

Although  the  model  described  above  accounts  for 
the  observed  correlations  between  Fos  expression  and 
phase  shifting,  it  is  important  to  emphasize  that  other 
possible  interpretations  of  these  data  are  equally  rea¬ 
sonable  and  that  thus  far  there  is  only  limited  direct 
evidence  for  a  causal  relationship  between  Fos  and 
phase  shifting  (Wollnik  et  al.,  1995;  Honrado  et  al., 
1996).  It  remains  to  be  discovered  which  genes  in 
which  SCN  cells  are  the  targets  of  transcriptional  regu¬ 
lation  by  c-fos  and  jun- B,  and  what  role  these  genes 
may  play  in  the  cellular  events  leading  to  resetting  the 
clock. 
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We  recorded  multiple  unit  neural  activity  [muitiunit  activity 
(MU A)]  from  inside  and  outside  of  the  suprachiasmatic  nucleus 
(SCN)  in  freely  moving  male  golden  hamsters  housed  in 
running-wheel  cages  under  both  light/dark  cycles  and  constant 
darkness.  The  circadian  period  of  MUA  in  the  SCN  matched  the 
period  of  locomotor  activity;  it  was  -24  hr  in  wild-type  and  20 
hr  in  homozygous  tau  mutant  hamsters.  The  peak  of  MUA  in  the 
SCN  always  occurred  in  the  middle  of  the  day  or,  in  constant 
darkness,  the  subjective  day.  There  were  circadian  rhythms  of 
MUA  outside  of  the  SCN  in  the  ventrolateral  thalamic  nucleus, 
the  caudate  putamen,  the  accumbens'  nucleus,  the  medial 
septum,  the  lateral  septum,  the  ventromedial  hypothalamic 
nucleus,  the  medial  preoptic  region,  and  the  stria  medullaris. 
These  rhythms  were  out-of-phase  with  the  electrical  rhythm  in 
the  SCN  but  in-phase  with  the  rhythm  of  locomotor  activity, 
peaking  during  the  night  or  subjective  night.  In  addition  to 
circadian  rhythms,  there  were  significant  ultradian  rhythms 
present;  one.  with  a  period  of  -80  min,  was  in  antiphase 


Circadian  locomotor  activity  rhythms  in  mammals  are  generated 
by  an  endogenous  pacemaker  located  in  the  suprachiasmatic 
nucleus  (SCN)  of  the  hypothalamus  (for  review,  see  Turek,  1985; 
Meijer  and  Rietveld.  1989;  Klein  et  al.,  1991).  Lesions  of  the  SCN 
cause  arrhv  thmicity  of  locomotor  activity  (Moore  and  Eichler, 
1972;  Stephan  and  Zucker,  1972;  Rusak  and  Zucker,  1979),  and 
transplants  of  fetal  SCN  tissue  restore  circadian  periodicities 
(Sawaki  et  al..  1984;  Lehman  et  al.,  1987;  Ralph  et  al.,  1990).  The 
SCN  exhibits  circadian  rhythms  in  several  in  vitro  preparations: 
the  acute  slice  (Green  and  Gillette,  1982;  Groos  and  Hendricks, 
1982;  Shibata  and  Moore,  1988),  slice  culture  (Bos  and  Mirmiran. 
1990;  Herzog  et  al.,  1997).  and  dispersed  cell  culture  (Welsh  et 
al..  1995:  Liu  et  al..  1997).  Both  slice  and  dispersed  cell  cultures  of 
SCN  also  display  circadian  rhythms  of  peptide  release  (Mu- 
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between  the  SCN  and  other  brain  areas,  and  another,  with  a 
period  of  -14  min,  was  in-phase  between  the  SCN  and  other 
brain  areas.  The  periods  of  these  ultradian  rhythms  were  not 
significantly  different  in  wild-type  and  tau  mutant  hamsters.  Of 
particular  interest  was  the  unique  phase  relationship  between 
the  MUA  of  the  bed  nucleus  of  the  stria  terminalis  (BNST)  and 
the  SCN;  in  these  two  areas  both  circadian  and  ultradian  com¬ 
ponents  were  always  in-phase.  This  suggests  that  the  BNST  is 
strongly  coupled  to  the  SCN  and  may  be  one  of  its  major  output 
pathways.  In  addition  to  circadian  and  ultradian  rhythms  of 
MUA,  neural  activity  both  within  and  outside  the  SCN  was 
acutely  affected  by  locomotor  activity.  Whenever  a  hamster  ran 
on  its  wheel,  MUA  in  the  SCN  and  the  BNST  was  suppressed, 
and  MUA  in  other  areas  was  enhanced. 

Key  words:  circadian;  ultradian;  suprachiasmatic  nucleus;  in 
vivo  recording;  hamster;  tau  mutant;  locomotor  activity;  bed 
nucleus  of  the  stria  terminalis;  MUA 


rakami  et  al.,  1991;  Watanabe  et  al..  1993;  Shinohara  et  al.,  1995). 
In  contrast  the  physiology  of  the  SCN  in  vivo  and  its  relationship 
to  circadian  behavior  in  the  intact  animal  have  received  little 
experimental  attention. 

To  understand  how  the  circadian  clock  controls  locomotor 
behavior,  we  need  to  understand  its  connections  to  the  motor 
control  system.  Although  output  pathways  from  the  SCN  circa¬ 
dian  pacemaker  are  not  completely  described,  the  motor  control 
system  in  mammals  is  relatively  well  characterized  (Wichmann  et 
al.,  1995;  Bergman  et  al.,  1998).  Because  there  are  no  known 
direct  neural  connections  between  the  SCN  and  motor  control 
areas  of  the  brain,  it  is  likely  that  either  humoral  factors  and/or 
“relay  nuclei”  serve  to  connect  the  SCN  with  the  motor  centers. 

Not  only  does  the  SCN  regulate  locomotor  activity  but  there  is 
reason  to  believe  that  locomotor  activity  feeds  back  on  the  SCN. 
The  period  of  the  free-running  rhythm  of  rats  housed  in  cages 
with  a  running  wheel  is  different  from  that  of  rats  housed  in  cages 
without  a  wheel  (Yamada  et  al.,  1988,  1990;  Shioiri  et  al.,  1990). 
Locking  of  the  running  wheel  changes  the  free-running  period  in 
mice  (Edgar  et  al.,  1991).  Access  to  running  wheels  induces  phase 
shifts  of  locomotor  activity  in  hamsters  (Mrosovsky,  1988;  Reebs 
and  Mrosovsky,  1989)  as  does  injection  of  triazolam,  which  in¬ 
creases  locomotor  activity  (van  Reeth  et  al.,  1987).  In  mice, 
forced  treadmill  running  induces  phase  shifts  of  circadian 
rhythms  and  is  able  to  entrain  them  (Marchant  and  Mistlberger, 
1996). 

The  apparent  complexity  of  the  relationship  between  the  SCN 
and  locomotor  centers,  almost  certainly  involving  reciprocal  in- 
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Figure  Examples  of  MUA  recorded  for  10  d  from  the  SCN  of  three 

4  "d  and  then  released  into  constant  darkness  (lighting  condition  in  ica  ,  •  tu  v_axis  for  this  portion  of  the  figure  shows  0-200 

Wheel-running  activity  is  plotted  at  the  bottom  of  A-C  as  t  e  num  er  o  revo  u  of  the  jeft  s^N  B  Recorded  from  the  ventromedial  portion  of 

from  the  ventromedial  portion  of  the  right  SCN  near  the  caudal  end  of  the  nucleus. 


teractions.  and  the  fact  that  direct  neuronal  interconnections 
appear  to  he  absent  provide  a  strong  rationale  for  exploring  the 
functional  relationships  in  vivo.  We  have  perfected  a  technique 
that  allows  u.s  to  record  neuronal  activity  of  several  brain  regions 
in  freely  moving  hamsters,  enabling  us  to  correlate  electrical 
activity  within  the  SCN  with  activity  in  other  brain  regions  and 
with  locomotor  activity.  We  have  used  this  technique  to  describe 
the  electrical  characteristics  of  the  SCN  in  vivo ,  the  differences 
between  the  tau  mutant  and  wild-type  hamsters,  the  relationship 
between  the  SCN  and  other  regions  of  the  brain,  and  the  effect  of 
the  animal's  locomotor  behavior  on  SCN  activity.  The  results 
provide  a  new  framework  for  understanding  the  regulation  of 
locomotor  behavior  by  the  circadian  timing  system. 

MATERIALS  AND  METHODS 

Animals.  Three-  to  five-month-old  golden  hamsters  (LVG  wild  type  from 
Charles  River  Laboratories.  Wilmington,  MA;  LVG  background  tau 
mutant  animals  from  our  colonv)  were  used  in  this  study.  Animals  were 
entrained  for  at  least  2  weeks  to  light/dark  cycles  (LDs)  (14:10  hr  LD  for 
wild  types;  1 1. 7:8.3  hr  LD  for  tau  mutants;  light  intensity  of  -dOO  lux  at 
cage  level).  We  monitored  wheel-running  activity  throughout  the  exper¬ 
iments  and  used  only  animals  that  showed  dear  locomotor  rhythmicity. 

Electrode  implantation.  We  implanted  one  or  two  bipolar  electrodes 
constructed  from  pairs  of  Teflon-coated  stainless  steel  wires  (bare  diam¬ 
eter  130  M-m;  A-M  Systems,  Everett,  WA;  tip  distance,  -150  /xm  for 
recording  from  the  SCN  and  200-300  for  other  brain  regions)  and  an 
uncoated  platinum-iridium  wire  (diameter,  130  ^m;  A-M  Systems)  used 
as  a  signal  ground  in  the  cortex.  Wires  were  connected  to  an  eight  pin  IC 
socket" wrapped  in  insulated  copper  tape.  Distances  between  any  two 
bipolar  electrodes  were  determined  according  to  the  recording  sites 

Electrode  implantation  was  performed  under  pentobarbital  anesthesia 
(90  mg/kg.  i.p.).  Animals  were  placed  in  a  stereotaxic  instrument 
with  the  nose  bar  set  at  -2  mm  (David  Kopf  Instruments,  Tujunga,  CA). 


Four  self-tapping  screws  (#0  x  1/8  inch;  Small  Parts,  Miami  Lakes,  FL) 
were  implanted  into  1  mm  holes  in  the  skull  made  with  a  dental  drill.  We 
used  different  stereotaxic  coordinates  for  wild-type  and  tau  mutant 
hamsters  because  the  shape  of  bregma  in  tau  mutants  is  different  and 
more  variable  than  is  that  in  wild  types.  Wild-type  SCN  coordinates  were 
0  7  mm  anterior  to  bregma,  0.2  mm  lateral  to  the  midsagittal,  and  8.0-8.2 
mm  below  the  dural  surface.  Tau  mutant  coordinates  were  10-1.2  mm 
anterior  to  bregma,  0.2  mm  lateral  to  the  midsagittal,  and  7.9-8.1  mm 
below  the  dural  surface.  The  electrode  was  secured  to  the  screws  and  the 
skull  with  dental  cement. 

Recording  procedure.  One  week  after  surgery,  each  hamster  was  trans¬ 
ferred  to  a  24  cm  (width)  x  21  cm  (length)  x  30  cm  (height)  cage  with 
a  running  wheel  21  cm  in  diameter  mounted  on  one  side  to  allow  the 
hamster  equipped  with  wires  access  to  the  wheel.  The  electrodes  were 
connected  to  head  stage  buffer  amplifiers  (J-FET  input  OP  Amp;  TL084) 
located  on  the  hamster’s  head.  Buffer  amplifiers  were  connected  to  a  12 
channel  slip  ring  (Airflyte  Electronics  Company,  Bayonne,  NJ)  that 
allowed  free  movement  for  the  animal.  The  wires  between  the  head  stage 
amplifiers  and  the  slip  ring  were  protected  by  a  stainless  steel  spring. 
Output  signals  were  processed  by  differential  input  integration  amplifiers 
(INA  101  AM;  Burr-Brown,  Tucson,  AZ;  gain,  x!0)  and  then  fed  into 
AC  amplifiers  (OP  Amp,  4558;  bandpass,  500  Hz  to  5  kHz;  gain,  10,000). 
Spikes  were  discriminated  by  amplitude  and  counted  in  1  min  bins  using 
a  computer-based  window  discrimination  system  (DAS-1801ST  AD 
board;  Keithley  Metrabyte,  Taunton,  MA).  Wheel  revolutions  were  re¬ 
corded  using  the  Data  Quest  system  (Data  Science  International,  St. 

Paul,  MN).  .  .  c  . 

Reduction  of  recording  noise .  In  recording  neural  activity  from  freely 
moving  animals,  the  biggest  problem  is  noise.  We  reduced  microphomc 
noise,  which  is  caused  by  mechanical  disturbances  such  as  wire  move¬ 
ments,  by  mounting  the  buffer  amplifiers  on  the  head  (effectively  decreas¬ 
ing  the  impedance).  We  directly  coupled  the  output  signals  to  the 
integration  amplifier  (without  any  capacitors  or  resistors)  to  provide  a 
high  common  mode  rejection  ratio  that  can  reduce  non-neuronal  signals 
such  as  muscle  potentials.  We  also  used  shielding  material  around  the 
electrode  and  its  vicinity  to  reduce  noise  from  the  animal  s  scratching. 
Because  we  could  not  entirely  remove  this  source  of  noise,  we  used  a 
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Figure  2.  Daily  and  circadian  rhythms  of  neural  activity  in 
several  regions  of  the  brain.  MUA  and  locomotor  activity 
are  plotted  in  6  min  bins  as  in  Figure  1.  Recordings  were 
made  from  the  following:  A,  right  side  of  the  ventrolateral 
thalamic  nucleus  (VLT)\  B ,  right  side  of  the  medial  septum 
(A/5);  C,  right  side  of  the  stria  medullaris  {sm):  D.  right  side 
of  the  LS;  E ,  the  optic  chiasm  (or).  Each  plot  represents  a 
different  wild-type  hamster  except  for  A  (same  animal  as  in 
E)  and  C  (same  animal  as  in  Fig.  1C).  Note  (most  clearly  in 
B)  that  the  peak  of  neural  activity  coincides  with  wheel¬ 
running  activity. 


highly  effective  low-cut  filter  (500  Hz).  Using  these  techniques,  we 
reduced  electrical  noise  generated  by  chewing  or  moving  to  undetectable 
levels.  Although  scratching  did  generate  detectable  electrical  noise,  this 
activity  was  rare  and  did  not  create  a  problem  in  the  analysis. 

Identification  of  recording  sites.  After  the  electrical  recordings,  each 
hamster  was  anesthetized  w  ith  halothanc,  the  head  amplifier  was  discon¬ 
nected,  and  a  small  positive  current  (50  /j.A;  10  sec)  was  passed  through 
the  recording  electrodes.  The  brain  was  removed  and  fixed  in  Zamboni's 
fixative  solution  for  a  few  days.  Frozen  sections  (40  ^.m  thick)  were 
stained  with  potassium  ferrocyanide  (5rr  potassium  ferrocyanide  in  lO'T 
HC1).  Blue  spots  of  deposited  iron  were  used  for  identification  of 
recording  sites. 

Data  analysis.  The  first  report  of  SCN  neuronal  recordings  from  freely 
moving  rodents  appeared  19  years  ago  and  played  a  pivotal  role  in 
identifying  the  SCN  as  the  central  mammalian  circadian  pacemaker 
(Inouve  and  Kawamura.  1979).  Since  then  little  use  has  been  made  of  this 
technique.  A  primary  reason  for  hesitation  in  the  use  of  in  vivo  recording 
techniques  is  that  changes  in  neuronal  oscillations  are  oftentimes  difficult 
to  identify  in  the  raw  data,  and  thus  robust  time  series  statistical  analysis 
is  required.  No  single  time  series  analysis  too!  can  be  applied  in  al!  cases. 
We  applied  a  new  method,  singular-spectrum  analysis  (SSA).  in  combi¬ 
nation  with  older  methods.  Periodicities  in  multiunit  activity  (MUA) 
recorded  in  vivo  were  determined  using  SSA  in  combination  with  the 
multitapcr  method  (MTM)  approach  to  the  fast  Fourier  transform 
(Thomson,  1982;  Vautard  et  al.,  1992). 

Singular-spectrum  analysis  or  SSA  is  a  linear,  nonparametric  method 
based  on  a  principal  component  analysis  in  the  vector  space  of  the  delay 
coordinates  for  a  times  series  (Eisner  and  Tsonis,  1996).  In  SSA,  a  single 
time  series  is  expanded  into  a  set  of  multivariate  time  series  of  length  A/, 
known  as  the  “window  length.”  M  determines  what  range  of  frequencies 


can  be  resolved  as  a  stationary  signal  in  the  calculated  principal  compo¬ 
nents.  The  principal  component  analysis  orders  the  expanded  time  series 
as  a  new  coordinate  system  with  most  information  along  the  first  coor¬ 
dinates.  The  principal  components  are  processes  of  length  N  -  M  +  1 
that  can  be  thought  of  as  weighted  moving  averages  of  the  time  series  in 
which  each  accounts  for  a  certain  percentage  of  the  total  variance.  SSA 
allows  optimal  detrending,  identification  of  the  noise  floor  in  spectral 
estimates,  and  identification  of  intermittent  oscillatory  components  in  the 
data.  In  practice,  we  found  that  we  could  with  reasonable  success  divide 
the  data  into  four  parts  (trend  variance,  circadian  variance,  ultradian 
components  variance,  and  noise  variance)  by  use  of  two  window  lengths 
on  the  MUA  data.  M  was  set  at  36  hr  ( M  —  360  for  6  min  bins)  for  the 
circadian  time  scale  and  5  hr  for  the  ultradian  time  scale.  SSA  cannot 
resolve  periods  longer  than  M  and  treats  them  as  trends.  If  M  is  much 
greater  than  the  average  lifetime  of  an  episode  of  oscillation,  SSA  cannot 
resolve  the  intermittent  oscillation. 

We  used  SSA  for  signal  reconstruction  from  the  noisy  MUA  data. 
Simple  noise  reduction  by  applying  a  fixed  low-pass  filter  to  the  data  is 
not  appropriate  when  the  spectrum  is  not  monotonic.  Because  steps  were 
taken  to  minimize  instrument  noise  and  a  low-cut  filter  was  used  (500 
Hz)  before  binning  the  impulses  in  the  MUA  data,  noise  is  represented 
here  mostly  by  random  impulses  from  populations  of  neurons  near  the 
electrode.  Optimal  filtering  of  signals  that  are  not  completely  stable 
requires  methods  such  as  Wiener  filtering  or  SSA.  Both  provide  optimal 
filters  in  a  least  squares  sense.  However  Weiner  filters  arbitrarily  require 
harmonic  functions  as  a  basis.  SSA,  in  contrast,  uses  data-determined 
general  functions  that  do  not  require  any  previous  hypotheses  about  the 
noise  variance.  Unlike  SSA,  the  Wiener  method  requires  smooth  and 
very  reliable  estimates  of  the  power  spectrum  that  are  impossible  to 
obtain  with  very  short  data  sets.  “Noise-free”  circadian  or  ultradian  time 
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series  were  therefore  reconstructed  from  the  SSA  filters.  By  this  method, 
oscillatory  signals  that  accounted  for  as  little  as  3 %  of  the  total  variance 
could  be  detected.  Monte  Carlo  simulations  assuming  either  white  noise 
or  correlated  noise  were  used  to  evaluate  the  background  noise  level  in 
the  SSA  spectra.  Bv  selecting  onlv  those  principal  components  associated 
with  the  circadian  variance  of  the  data,  we  could  derive  an  optimal  (in 
the  least  squares  sense)  -  noiseless”  reconstruction  of  the  circadian  wave¬ 


form  in  the  data.  Each  waveform  in  the  reconstructed  time  series  is 
essentially  a  local  fit  allowing  us  to  calculate  a  mean  period  from  each 
animal's  data  set.  This  local  noise-free  fit  also  allowed  us  to  determine 
the  phase  relationship  as  a  function  of  time  between  two  different 
locations  in  the  brain.  Periods  of  the  ultradian  rhythms  were  determined 
with  one  or  more  of  the  following:  SSA-MTM,  power  spectrum  estimates 
[e.g.,  maximum  entropy  method  (M  EM)],  or  visual  inspection  of  the  raw 


WT#1 


Time  in  DD  (h) 


FiS„r,  *  Ultradian  activity  record,  from  the  SC N and  LS  ol  ■ r.M-W  hamster^.  SmSHm  Xd 

correlation  of  neural  activity  with  wheel-running  activity  in  the  SCN  and  the  pos.tive  correlation  in  the  LS. 
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Figure  9  Mathematical  analysts  of  ultradian  rhythms  recorded  from  the  SCN  and  LS  (from  Fig.  3A.B).  A.  Data  reconstructed  by  SSA  of  the  80  min 
ultradian  rhythms  in  the  SC  N  (solid  line)  and  the  LS  (dotted  line)  (ox  the  first  72  hr  of  the  record.  B,  Phase  angle  difference  of  the  80  nun  ultradian  rhythms 
between  the  SCN  and  the  LS  plotted  for  144  hr  (see  text  and  Fig.  4  for  details). 


data.  MEM  was  performed  on  the  data  after  the  trend  and  circadian 
variances  were  removed  (determined  by  SSA).  The  order  of  the  MEM 
(Af.  number  of  poles  in  the  autoregressive  model)  was  kept  much  lower 
than  the  number  of  data  points  and  varied  over  the  range  to  monitor 
stability  of  peaks  of  interest  in  the  spectrum  (Af,  20-60;  N  >  1600). 
MEM  is  fully  consistent  with  SSA  so  that  the  additive  property  of  the 
spectra  is  conserved  exactly  (Vautard  et  al.,  1992).  Periods  of  very  short 
ultradian  rhythms  were  determined  by  visual  inspection;  areas  of  high- 
frequency  ultradian  activity  were  measured  for  the  average  period,  and 
measurements  from  six  areas  per  animal  were  averaged  to  obtain  the 
period.  The  periods  of  the  wheel-running  activity  were  obtained  by  using 


a  standard  chi-square  periodogram  for  comparison  with  the  other  meth¬ 
ods  (see  Table  1). 

Phase  analysis.  We  used  the  discrete  Hilbert  transform  of  the  SSA- 
reconstructed  waveforms  to  estimate  the  instantaneous  phase  difference 
between  the  oscillations  inside  and  outside  the  SCN.  The  Hilbert  transform 
is  the  imaginary  part  of  the  analytical  signal  of  a  real  time  series  (Bendat 
and  Piersol,  1986).  An  estimate  of  the  imaginary  part  from  the  real  part  of 
the  time  series  can  be  calculated  by  assuming  the  signal  is  causal  (a 
sequential  time  series)  and  by  using  the  inverse  of  a  one-sided  Fourier 
transform  (fast  Fourier  transform  coefficients  that  correspond  to  negative 
frequencies  are  set  to  zero  before  doing  the  inverse  transform).  The 
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Ultradian  rhvthms  in  neural  activity  records  from  the  SCN  and  CP  in  a  tau  mutant  hamster.  .4,  Expanded  plot  of  24  hr  neural  activity  in 
SCN  (from  Fig.  5,4,  double-headed  arrow).  The  top  bar  shows  wheel-running  activity,  plotted  as  a  black  bar  whenever  the  wheel 
B.  Expanded  plot  of  the  concurrent  record  from  the  CP  (from  Fig.  5B ,  double-headed  arrow).  Note  the  positive  (CP) 


Figure  10. 

T  min  bins  from  the 
revolved  more  than  once  in  6  min. 

and  negative  (.S'C.Y)  correlations  with  wheel-running  behavior. 

resultant  time  series  is  phase  shifted  90°  from  the  original  time  series.  This 
process  is  sensitive  to  noise  and  requires  a  filtered  input  such  as  that 
provided  bv  SSA.  The  temporal  change  in  the  relative  phase  between  the 
oscillations  recorded  in  different  brain  regions  was  examined  using  this 
method.  Circular  statistics  were  used  to  compute  the  mean  phase  difference 
for  the  circadian  rhythms  in  different  brain  regions  (Fisher,  1995). 

Experimental  protocols.  To  observe  daily  and  circadian  changes  of  neural 
activity  in  and  outside  of  the  SCN.  we  recorded  MUA  from  20  wild-type 
hamsters  implanted  with  two  electrodes.  One  electrode  w'as  directed  into 
the  SCN.  and  the  other  was  aimed  at  one  of  several  regions  outside  of  the 
SCN.  Animals  were  recorded  tor  4  d  in  a  light  dark  cycle  (14: 10  hr  LD;  light 
intensity  ot  -300  lu\  at  cage  level),  followed  by  6  d  in  constant  darkness 
(DD).  To  evaluate  the  effects  ot  the  tau  mutation  on  the  rhythmic  proper¬ 
ties.  we  recorded  MI  A  from  111  wild-type  and  14  tau  mutant  hamsters 
implanted  with  two  electrodes,  one  aimed  at  the  SCN  and  the  other  at  one 
of  the  following  sites:  the  Literal  septum  (LS).  the  caudate  putamen  (CP), 
the  ventrolateral  region  of  the  thalamus,  and  the  bed  nucleus  of  the  stria 


terminalis  (BNST).  These  hamsters  were  placed  into  DD  at  the  time  of 
their  normal  lights  off,  and  MUA  was  recorded  for  7  d.  We  only  used  data 
from  implantations  in  which  both  wires  of  each  electrode  were  located  in 
the  same  nucleus.  We  sometimes  observed  one  wire  located  in  the  SCN, 
while  the  other  was  in  the  surrounding  area;  in  those  cases  the  data  w'ere 
discarded. 

RESULTS 

Circadian  rhythms  both  within  and  outside  the  SCN 

MUA  recorded  in  the  SCN  showed  clear  daily  (LD)  and  circadian 
(DD)  rhythms  (Fig.  1).  Peak  impulse  frequency  of  these  rhythms 
always  occurred  at  the  middle  of  the  day  (LD)  or  subjective  day 
(DD),  in  antiphase  with  the  hamster’s  nocturnal  wheel-running 
activity.  The  amplitude  of  the  neuronal  activity  rhythms  varied 
from  experiment  to  experiment  and  was  not  strongly  correlated 


Table  1.  Summary  of  circadian  and  ultradian  periods 


Wheel-running  activity 
(hr) 


WT#  1 

24.3 

wt#: 

24.2 

WT#3 

23.9 

WT#5 

24.2 

Mean  =  SE 

' 

24.15  =  0.09 

ss#: 

19.8 

SS#4 

20.1 

SS#5 

20.0 

SS#7 

19.8 

Mean  *  SE 

19.93  ±  0.08 

Two-tailed  t  test 

p  <  0.0005 

Neural  activity  in  the  SCN 


Circadian 

(hr) 

Ultradian  1 
(min) 

Ultradian  ! 
(min) 

24.5 

84 

14 

24.1 

75 

14 

23.8 

87 

14 

24.1 

76 

15 

24.13  ±  0.14 

80.5  ±  3.0 

14.3  ±  0.3 

19.7 

75 

15 

20.3 

78 

15 

20.1 

99 

12 

19.7 

71 

15 

19.95  ±  0.15 

80.8  ±  6.3 

14.3  ±  0.8 

p  <  0.0005 

p  =0.97 

p  =  1-00 

Seven  davs  of  MUA  and  wheel-running  activity  were  collected  simultaneously  in  constant  darkness  (WT,  wild  type;  SS,  homozygous  tau  mutants)  Periods  of  wheel-running 
rhythms  were  obtained  by  chi-square  periodogram.  Periods  of  circadian  and  ultradian  l  rhythms  of  MUA  were  determined  by  SSA.  Periods  of  ultradian  2  rhythms  were 
determined  by  visual  inspection.  (See  text  for  details  of  analysis  methods.) 
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SS#2 


Figure  11.  Mathematical  analysis fwlhefim  72hr“  a^gle  difference  of  the  80  min 

There  was  no  consistent  relationship  between  overall  levels  of 
locomotor  activity  and  MUA  activity.  For  example,  the  wheel¬ 
running  activity  of  animals  shown  in  Figure  l,  A  and  B,  was  quite 
variable  (on  some  nights  the  hamsters  ran  much  more  than  on 
others),  yet  the  amplitude  of  the  MUA  rhythms  remained  rela- 
lively  constant.  Furthermore,  in  all  three  examples  in  Figure  1, 
the  amplitude  of  MUA  oscillation  in  the  SCN  did  not  change 
when  the  animals  were  released  from  LD  into  DD.  Even  in  the 
case  of  the  animal  whose  data  are  shown  in  Figure  IA  m  which 
the  electrode  was  located  in  the  lateral  portion  of  the  SCN  (which 


with  the  recording  site,  although  recordings  obtained  in  the 
ventral  portion  of  the  SCN  tended  to  have  a  higher  amplitude 
than  did  those  from  the  dorsal  portion.  Although  the  recording 
sites  were  nearly  identical  for  the  animals  whose  records  are 
shown  in  Figure  1,  B  and  C,  the  MUA  displayed  in  Figure  IB  has 
higher  average  amplitude  oscillations  than  does  that  in  Figure  1C. 
SSA  detected  a  significant  circadian  (or  daily)  oscillation  with  a 
centered  amplitude  of  -10.000  spikes  per  6  min  interval  from  the 
data  plotted  in  Figure  IB.  The  centered  amplitude  from  data 
shown  in  Figure  1C  was  approximately  one-fourth  as  large. 
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Figure  12.  Neural  activity  records  obtained the^ighTanterom^dL1  portion  of  the  BNST.  A, 
located  in  the  right  ventromedial  jwrtion  of  marked  by  the  double-headed  arrows  are 


constant  darkness.  Electrodes  were 
MUA  in  the  SCN  in  6  min  bins.  B, 
presented  below  (see  Fig.  14).  Note 


(^^and  for^the1  BNST^B^c'lnstarUaneous  phase  angle  difference  between  the 

circadian  rhythms  shown  in  A  and  B  (see  text  and  Fig.  4  for  details). 
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Figure  14.  Uitradian  rhvthms  of  neural  activity  recorded  from  the  SCN  and  the  BNST  in  a  wild-type  hamster.  A,  Expanded  plot  of  24  hr  neural  activity 
record  in  1  min  bins  from  the  SCN  (from  Fig.  12,4,  double -headed  arrow).  The  top  bar  shows  wheel-running  activity,  plotted  as  a  black  bar  whenever  the 
wheel  revolved  more  than  once  in  6  min.  5,  Expanded  plot  of  24  hr  neural  activity  record  in  1  min  bins  from  the  BNST  (from  Fig.  125,  double-headed 
arrow).  Note  that  both  the  80  and  14  min  uitradian  rhythms  are  usually  in-phase. 


receives  the  retinal  projection),  light  did  not  have  an  acute  effect 
on  the  MUA. 

Most  regions  outside  of  the  SCN  also  showed  clear  daily  or 
circadian  rhythms  (Fig.  2).  The  ventrolateral  thalamic  nucleus 
(n  =  3),  the  CP  in  =  3),  the  accumbens  nucleus  (/j  =  l),  the 
medial  septum  (n  =  1),  the  LS  (n  =  6),  the  ventromedial 
hypothalamic  nucleus  (n  =  2),  the  medial  preoptic  region  (n  = 
1).  and  the  stria  medullaris  (n  =  3)  all  exhibited  circadian 
rhythms  with  peak  impulse  activity  occurring  at  night  or  sub¬ 
jective  night  that,  unlike  the  electrical  activity  rhythm  in  the 
SCN.  was  in-phase  with  locomotor  activity.  SSA  of  activity  in 
the  optic  chiasm  (//  =  2)  revealed  no  circadian  components,  but 
uitradian  rhythms  were  present  and  were  stronger  in  LD  than 
in  DD. 

There  were  clear  correlations  between  electrical  activity  in 
the  SCN  and  in  other  brain  regions.  The  data  in  Figure  3  were 
recorded  from  a  w  ild-type  hamster  in  DD.  MUA  from  the  SCN 
(Fig.  3.-1 )  and  the  LS  (Fig.  35)  and  wheel-running  activity  (Fig. 
3C)  were  collected  simultaneously;  all  three  data  sets  contain 
significant  circadian  components.  Figure  4  presents  recon¬ 
structed  waveforms  of  the  circadian  components  using  the  SSA 
method.  Circadian  rhythms  in  the  SCN  (Fig.  4 A)  and  LS  (Fig. 
45)  were  tightly  locked  in  an  antiphase  relationship  (Fig.  4C). 
Their  periods  were  matched  to  that  of  the  wheel-running 
activity  rhythm  that  was  in  antiphase  to  the  SCN.  The  data  in 
Figure  5  were  recorded  from  a  tau  mutant  hamster.  MUA  was 
recorded  from  the  SCN  (Fig.  5/1)  and  CP  (Fig.  55).  Both 
regions  showed  significant  circadian  rhythms  that  matched  the 
20  hr  period  of  the  wheel-running  rhythm  (Fig.  5C).  The  phase 
angle  difference  between  these  two  brain  regions  was  also 
stable  in  antiphase  (Fig.  6;  sometimes,  as  in  this  example, 
stable  phase  relationships  did  not  develop  the  first  cycle  or  two 
of  constant  darkness).  Other  examples  (SCN  and  LS  in  a  wild 
type,  n  =  1;  SCN  and  ventrolateral  thalamus  in  a  wild  type, 


n  -  1;  and  SCN  and  ventrolateral  thalamus  in  a  tau  mutant, 
n  =  1)  showed  similar  phase  relationships. 

Uitradian  rhythms 

MUA  in  the  SCN  and  in  other  brain  regions  displayed  uitradian 
components  in  addition  to  circadian  ones.  We  succeeded  in 
identifying  a  significant  circadian  component  and  two  uitradian 
components  from  the  time  series  data  in  all  brain  regions  ana¬ 
lyzed.  Figure  7  shows  an  example  of  MUA  recorded  from  the 
SCN  in  a  tau  mutant  hamster  in  DD.  As  expected,  an  -20  hr 
circadian  component  in  antiphase  to  the  wheel-running  activity 
rhythm  was  present.  The  record  also  shows  a  clear  uitradian 
component  with  a  period  of  -80  min.  Furthermore,  in  expanded 
plots  we  could  detect  a  higher  frequency  rhythm  of  —14  min  that 
was  most  clearly  expressed  around  the  peak  of  the  80  min 
oscillation. 

Expanded  plots  of  24  hr  segments  recorded  from  the  SCN  (Fig. 
3/1)  and  the  LS  (Fig.  35)  of  a  wild-type  hamster  are  shown  in 
Figure  8,  A  and  5,  respectively.  These  data  contain  two  uitradian 
components,  -80  and  -14  min.  The  two  components  were  iden¬ 
tified  by  spectral  analysis  and  visual  inspection,  respectively.  The 
time  domain  SSA  reconstructions  of  the  —80  min  uitradian  oscil¬ 
lations  in  the  SCN  and  the  LS  were  in-phase  with  each  other  on 
the  first  day  after  transfer  from  LD  to  DD;  however  by  the  second 
day  in  DD  they  had  shifted  to  an  antiphase  relationship  (Fig. 
9/1,5).  The  —14  min  uitradian  rhythms  appeared  to  be  in-phase 
between  the  SCN  and  the  LS  (Fig.  8).  We  obtained  similar  results 
from  two  other  wild-type  hamsters  (recorded  from  the  SCN  and 
the  ventrolateral  thalamus  or  LS). 

The  uitradian  components  recorded  from  tau  mutant  hamsters 
also  had  periods  of  -80  and  -14  min.  Time  series  analysis  and 
visual  inspection  of  simultaneous  recordings  from  the  SCN  and 
the  CP  in  a  tau  mutant  hamster  revealed  these  two  uitradian 
components  (Fig.  10).  In  this  data  set  the  phase  relationships  of 
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the  80  min  oscillations  in  the  two  areas  were  often  in  antiphase, 
but  their  phase  relationship  was  much  less  stable  than  in  similar 
recordings  from  wild  types  (Fig.  11).  As  in  wild  types  the  ~14 
min  rhythms  always  appeared  to  be  in-phase  (Fig.  10).  A  second 
tau  mutant  hamster  in  which  electrical  recordings  were  obtained 
from  the  SCN  and  the  ventrolateral  thalamus  exhibited  simi  ar 
properties. 


Effects  of  the  tau  mutation  on  the  circadian  and 
ultradian  periods  of  MUA  . 

The  period  of  the  rhythm  of  wheel-running  activity  in  wild-type 
hamsters  was  ~24  hr  and  closely  matched  the  circadian  compo¬ 
nent  of  the  electrical  activity  in  the  SCN.  To  assess  the  effect  o 
the  tau  mutation  on  the  rhythmic  properties  in  the  SCN  we 
obtained  four  1-week-long  recordings  for  analysis  of  the  period  o 
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the  rhythm  in  DD.  The  ~20  hr  period  of  wheel-running  rhythms 
in  tau  mutants  closely  matched  the  circadian  MUA  rhythms  in  the 
SCN.  There  were  no  significant  differences  between  wild  type 
and  tau  mutants  in  the  periods  of  the  ultradian  components  (80 
and  14  min).  The  variation  in  frequency  among  animals  within 
each  genotype  was  as  great  as  that  between  genotypes  (Table  1). 

Effect  of  wheel  running  on  MUA 

Locomotor  activity  affects  the  expression  of  MUA  in  the  SCN 
and  elsewhere  in  the  brain.  MUA  in  the  SCN  was  decreased 
during  wheel  running,  whereas  MUA  in  other  areas  was  enhanced 
(Figs  *1-3,  5,  7,  8.  10).  These  changes  in  the  levels  of  MUA  were 
precisely  correlated  with  locomotor  activity  and  were  more  pro¬ 
nounced  the  more  vigorously  the  animals  ran  in  their  wheels. 

Phase  of  MUA  in  the  BNST 

In  contrast  to  MUA  from  all  other  recording  sites  outside  of  the 
SCN.  MUA  recorded  from  the  anteromedial  part  of  the  BNST 
exhibited  oscillations  (circadian  and  ultradian)  in-phase  with  the 
SCN  (Figs.  12.  13,  14,  15).  As  seen  in  the  SCN  and  again  in 
contrast  to  that  seen  at  other  sites,  neural  activity  in  the  BNST 
was  suppressed  during  wheel  running  (Figs.  12,  14).  We  also 
observed  similar  phase  relationships  in  one  recording  from  the 
posteromedial  part  of  the  BNST  of  a  wild-type  animal  and  from 
the  anteromedial  part  of  the  BNST  in  one  tau  mutant. 

DISCUSSION 

Circadian  rhythms  in  electrical  activity 

Inouye  and  Kawamura  (1979)  reported  that  the  MUA  in  the  SCN 
recorded  trom  treelv  moving  rats  displayed  a  circadian  rhythm 
with  peak  impulse  frequency  occurring  during  the  day.  They  also 
recorded  MUA  circadian  rhythms  from  regions  outside  the  SCN 
that  exhibited  peak  activity  at  night,  out-of-phase  with  the  SCN 
rhythm  and  in-phase  with  the  animal’s  locomotor  activity.  Be¬ 
cause  MUA  in  the  SCN  remained  rhythmic  after  the  nucleus  was 
isolated  (within  a  hypothalamic  island),  while  the  rhythms  out¬ 
side  of  the  SCN  disappeared,  they  concluded  that  the  circadian 
pacemaker  was  in  the  SCN  and  that  rhythms  outside  the  SCN 
were  beinu  driven  (in  antiphase)  by  the  SCN  (Inouye  and 
Kawamura'.  19X2).  Later.  Yamazaki  et  al.  (1994)  showed  that  the 
circadian  rhythm  of  ATP  content  in  the  SCN  of  rats  is  also  in 
antiphase  with  this  rhythm  in  the  anterior  hypothalamus  and  CP. 
Our  results  with  hamsters  are  in  good  general  agreement  with  the 
results  in  rats.  We  find  that  the  circadian  rhythms  of  MUA  within 
the  SCN  are  out-of-phase  with  those  recorded  from  brain  areas 
outside  the  SCN,  with  the  important  exception  of  the  BNST,  an 
area  from  which  Inouye  and  Kawamura  did  not  record. 

We  found  that  the  period  of  the  circadian  rhythm  in  MUA  both 
within  and  outside  of  the  SCN  was  affected  as  expected  by  the  tau 
mutation  and  matched  the  period  of  wheel-running  behavior  both 
in  wild-type  hamsters  (~24  hr)  and  in  tau  mutants  (~20  hr). 
Davies  and  Mason  (1994)  reported  a  circadian  electrical  activity 
rhythm  with  a  period  of  —20  hr  from  tau  mutant  SCN  in  acute 
slice  preparation  that  was  monitored  for  36  hr.  Liu  et  al.  (1997) 
reported  that  circadian  periodicities  of  electrical  activity  in  SCN 
cells  dispersed  in  cell  culture  from  homozygous,  heterozygous, 
and  wild-type  hamsters  exhibit  free-running  periods  that  corre¬ 
late  with  their  genotype.  Our  results  obtained  in  vivo  confirm  that 
the  altered  period  of  neuronal  activity  is  expressed  within  the 
SCN  of  intact  tau  mutant  hamsters.  The  correspondence  of  the  in 
vivo  with  the  in  vitro  results  provides  important  confirmation  of 


the  usefulness  of  the  in  vitro  model.  We  know  that  at  least  the 
normal  period  of  the  intact  SCN  is  conserved  in  the  dish. 

The  genetically  determined  period  of  the  SCN  is  also  con¬ 
served  in  transplantation  experiments  in  which  the  locomotor 
rhythms  restored  by  transplanted  SCN  tissue  have  the  periods 
determined  by  the  genotype  of  the  donor  not  the  host  (Ralph  et 
al.,  1990;  Vogelbaum  and  Menaker,  1992;  Silver  et  al.,  1996). 
However,  the  data  of  Silver  et  al.  (1996)  suggest  strongly  that  in 
the  transplant  experiments  the  link  between  the  transplanted 
tissue  and  the  host’s  locomotor  system  is  hormonal  not  neural. 
Therefore,  although  the  several  studies  of  SCN  rhythmicity  using 
in  vitro,  transplant,  and  in  vivo  methods  are  consistent,  they  leave 
open  questions  about  the  functional  connections  of  the  SCN  with 
the  motor  control  system.  In  particular,  it  will  be  important  to 
determine  the  roles  of  both  neural  and  humoral  SCN  outputs  in 
controlling  locomotor  behavior  in  intact  animals. 

Ultradian  rhythms 

MUA  from  the  SCN  exhibited  two  different,  clearly  discernable 
ultradian  components  after  several  days  in  constant  darkness.  The 
80  min  ultradian  rhythm  was  out-of-phase  with  other  brain  re¬ 
gions  (except  the  BNST  in  which  all  rhythms  were  always  in-phase 
with  those  in  the  SCN).  However,  immediately  after  transfer 
from  LD  to  DD,  the  80  min  rhythms  within  and  outside  the  SCN 
were  temporarily  in-phase.  In  the  tau  mutant,  these  rhythms 
sometimes  ran  independently.  This  suggests  that  80  min  ultradian 
pacemakers  may  be  located  both  in  and  outside  of  the  SCN  and 
that  the  light/dark  cycle  may  influence  their  phase  relationships. 
A  similar  ultradian  rhythm  can  be  observed  in  figures  published 
by  Kawamura  and  Inouye  (1979;  their  Fig.  5)  and  Inouye  and 
Kawamura  (1982;  their  Fig.  7),  which  show  data  from  MUA 
recorded  in  rat  within  a  hypothalamic  island  containing  the  SCN 
as  well  as  outsi'de  of  the  island.  That  observation  supports  the  idea 
that  ultradian  periodicities  are  generated  from  sites  both  within 
and  outside  the  SCN.  Recently,  Meijer  et  al.  (1997)  found  signif¬ 
icant  ultradian  periods  of  4  hr  and  170.  130,  and  100  min  in  rat 
SCN  in  vivo  and  a  significant  period  of  3.5  hr  in  rat  SCN  in  vitro. 
MUA  from  optic  chiasm  in  hamster  and  rat  (Inouye  and 
Kawamura,  1979;  Omata  and  Kawamura,  1988)  has  been  re¬ 
ported  to  contain  circadian  components.  However,  in  our  record¬ 
ings  we  observed  only  ultradian  rhythms  from  this  site. 

The  source  and  mechanism  for  the  generation  of  ultradian 
periodicities  within  the  brain  are  obscure,  and  there  are  many 
possible  candidates  [e.g.,  in  the  posterior  hypothalamic  area  of 
the  conscious  rat,  the  release  rates  of  the  catecholamines  and 
histamine  fluctuate  with  the  following  periods:  histamine,  83  min; 
dopamine  and  noradrenaline,  92  min;  and  adrenaline,  99  min 
(Dietl  et  al.,  1992;  Prast  et  al.,  1992;  Grass  et  al.,  1996)]. 

Very  short  (-14  min)  ultradian  rhythms  were  present  in  all 
brain  regions  from  which  we  recorded.  Although  not  observable 
by  our  spectral  analysis  methods  because  of  limitations  in  sam¬ 
pling  frequency  resolution  and  signal  frequency  variability,  this 
periodicity  was  easily  detectable  visually  in  displays  of  the  record¬ 
ings.  The  rhythm  was  in-phase  among  all  brain  regions.  The  -14 
min  period  is  significantly  longer  than  an  ultradian  periodicity 
reported  previously  by  Miller  and  Fuller  (1992),  who  observed 
rhythms  with  periods  of  -120  sec  that  were  lengthened  by  retinal 
illumination  in  urethane-anesthetized  rats. 

We  failed  to  find  any  differences  in  the  periods  of  ultradian 
rhythms  of  tau  mutant  and  wild-type  hamsters,  even  though  the 
circadian  periodicity  of  the  tau  mutant  is  shortened  by  >15 
percent.  This  suggests  that  the  tau  mutation  may  affect  circadian 
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but  not  ultradian  rhythms  of  MUA.  This  is  in  contrast  to  its  effect 
on  luteinizing  hormone  pulsatility  (Loudon  et  al  1994)  and  to 
the  effect  of  mutations  in  the  per  gene  in  Drosophila  that  alter  the 
periods  of  both  circadian  rhythms  and  high-frequency  neural 
rhvthms  involved  in  the  control  of  the  courtship  song  (Kynacou, 
1990)  However,  our  data  do  suggest  that  ultradian  rhythms  of 
MUA  in  tau  mutants  may  have  more  variable  periods  and  phase 
relationships  than  do  those  in  wild  types,  although  we  are  pres¬ 
ently  unable  to  quantify  this  difference. 

Bed  nucleus  rhythmicity 

Unexpectedly,  and  uniquely  among  all  sites  in  which  record¬ 
ings  have  been  made,  the  circadian  rhythm  of  electrical  activity 
in  the  BNST  is  in-phase  with  the  SCN  electrical  activity 
rhythm.  In  a  previous  study  Inouye  (1983)  recorded  MUA  from 
seven  regions  of  the  rat  brain  (not  including  the  BNST),  but 
only  the  SCN  showed  a  circadian  rhythm  that  peaked  during 
the'dav:  rhvthms  from  all  the  other  sites  peaked  at  night.  In  our 
recordings!  ultradian  rhythms  in  the  BNST  were  also  in-phase 
with  those  in  the  SCN.  Cross-spectral  estimates  suggest  that 
the  BNST  and  SCN  are  tightly  coupled  to  each  other  (data  not 

shown).  ,  u 

There  is  anatomical  evidence  that  SCN  efferent  pathways  reach 
the  BNST  (Kalsbeek  et  al..  1993;  Morin  et  al.,  1994).  However,  it 
does  not  appear  that  the  BNST  projects  to  the  SCN  (Numan  and 
Numan.  1996).  The  BNST  may  therefore  be  part  of  an  output 
pathway  from  the  SCN  to  locomotor  centers.  This  would  be 
completely  consistent  with  our  data.  The  BNST  is  a  complex 
nucleus  with  several  subdivisions;  it  has  been  reported  to  be 
involved  in  the  control  of  photoperiod  measurement  (Raitiere  et 
al.,  1995),  maternal  behavior  (Numan  and  Numan,  1996),  an 
mating  behavior  (Wood  and  Newman,  1993),  all  of  which  may  be 
influenced  bv  the  circadian  system.  Our  preliminary  data  from 
lesion  studies  suggest  that  the  BNST  may  be  involved  in  control¬ 
ling  the  level  and  pattern  of  locomotor  activity  (Yamazaki  et  al., 
1997). 

Effects  of  locomotor  activity  on  SCN  electrical  activity 

Our  data  demonstrate  that  an  animal’s  movements  affect  electri¬ 
cal  activity  in  several  brain  areas.  Wheel-running  activity  acutely 
decreased  SCN  (and  BNST)  neural  activity  and  enhanced  neural 
activity  outside  the  SCN  (see  also  Meijer  et  al.,  1997).  Although 
we  were  not  able  to  determine  the  causal  relationship  involved 
(j.e..  does  increased  activity  in  extra-SCN  regions  lead  to  inhibi¬ 
tion  of  SCN  activity  or  vice  versa?),  we  believe  that  this  finding 
has  great  potential  significance.  It  is  intriguing  to  speculate  that 
reductions  in  SCN  neural  activity  associated  with  locomotor 
activity  mav  underlie  the  "nonphotic"  phase  shifts  of  circadian 
rhythms  that  are  produced  by  vigorous  wheel  running  (Reebs  and 
Mrosovskv,  1989). 

The  work  reported  here  represents  a  level  of  analysis  of  mam¬ 
malian  circadian  organization  that  has  been  almost  completely 
neglected  since  it  was  pioneered  by  Inouye  and  Kawamura  20 
years  ago.  Technical  and  conceptual  advances  over  the  past  two 
decades  have  made  it  more  tractable,  although  it  is  still  labor 
intensive.  In  spite  of  the  inherent  difficulty  of  obtaining  it,  knowl¬ 
edge  of  the  dynamics  of  electrical  activity  in  the  several  brain 
regions  controlling  circadian  behavior  in  intact,  unanesthetized 
animals  is  likely  to  be  essential  for  an  understanding  of  this 
important  regulatory  system. 
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Abstract  Circadian  pacemakers  in  many  animals  are  compound.  In  rodents,  a 
two-oscillator  model  of  the  pacemaker  comprising  an  evening  (E)  and  a  morning 
(M)  oscillator  has  been  proposed  based  on  the  phenomenon  of  "splitting"  and 
bimodal  activity  peaks.  The  authors  describe  computer  simulations  of  the  pace¬ 
maker  in  tau  mutant  hamsters  viewed  as  a  system  of  mutually  coupled  E  and  M 
oscillators.  These  mutant  animals  exhibit  normal  type  1  PRCs  when  released  into 
DD  but  make  a  transition  to  a  type  0  PRC  when  held  for  many  weeks  in  DD.  The 
two-oscillator  model  describes  particularly  well  some  recent  behavioral  experi¬ 
ments  on  these  hamsters.  The  authors  sought  to  determine  the  relationships 
between  oscillator  amplitude,  period,  PRC,  and  activity  duration  through  com¬ 
puter  simulations.  Two  complementary  approaches  proved  useful  for  analyzing 
weakly  coupled  oscillator  systems.  The  authors  adopted  a  "distinct  oscillators" 
view  when  considering  the  component  E  and  M  oscillators  and  a  "system"  view 
when  considering  the  system  as  a  whole.  For  strongly  coupled  systems,  only  the 
system  view  is  appropriate.  The  simulations  lead  the  authors  to  two  primary 
conjectures:  (1)  the  total  amplitude  of  the  pacemaker  system  in  tau  mutant  ham¬ 
sters  is  less  than  in  the  wild-type  animals,  and  (2)  the  coupling  between  the  unit  E 
and  M  oscillators  is  weakened  during  continuous  exposure  of  hamsters  to  DD. 
As  coupling  strength  decreases,  activity  duration  (a)  increases  due  to  a  greater 
phase  difference  between  E  and  M.  At  the  same  time,  the  total  amplitude  of  the 
system  decreases,  causing  an  increase  in  observable  PRC  amplitudes.  Reduced 
coupling  also  increases  the  relative  autonomy  of  the  unit  oscillators.  The  rela¬ 
tively  autonomousphaseshiftsof  E  and  M  oscillators  can  account  forboth  imme¬ 
diate  compression  and  expansion  of  activity  bands  in  tau  mutant  and  wild-type 
hamsters  subjected  to  light  pulses. 


Key  words  circadian  rhythm,  computer  simulation,  limit  cycles,  oscillators,  coupling, 
period  mutants 


Two  properties  of  the  circadian  pacemaker  play 
central  roles  in  the  mechanism  of  its  entrainment  by 
light-dark  cycles:  its  free-running  period  (x)  and  its 
phase  response  curve  (PRC).  The  PRC  shape  and 
amplitude  define  the  phase  relationship  attained  by 
an  oscillator  of  period  xfora  given  lightcycle  and  also 


determine  its  limits  of  entrainment  (Daan  and 
Pittendrigh,  1976).  Organisms  with  x  altered  by 
genetic  mutations  offer  a  good  opportunity  to 
approach  an  explanation  of  the  mechanistic  intercon¬ 
nections  between  x  and  the  PRC  (Menaker,  1992; 
Menaker  and  Takahashi,  1995). 
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Comparative  PRCs  for  period  mutant  organisms 
were  performed  first  for  eclosion  rhythms  (Konopka, 

1979)  and  later  for  activity  rhythms  (Saunders  et  al., 
1994)  in  Drosophila  and  in  Neurospora  (Dharmananda, 

1980) .  In  Drosophila,  the  per  short  mutation  led  to  a 
change  in  the  PRC  from  a  low-amplitude  type  1  to  a 
high-amplitude  type  0.  In  Neurospora,  different  period 
mutants  presented  a  correlation  between  period 
decreases  and  phase  response  increases,  with  the 
short  period  mutants  presenting  type  0  PRCs.  Com¬ 
parative  studies  of  PRCs  in  wild-type  (t  =  24  h)  and 
homozygous  tau  mutant  (x  =  20  h)  hamsters  were  per¬ 
formed  by  Menaker  and  his  coworkers  (Shimomura 
and  Menaker,  1994;  Menaker  et  al.,  1994;  Menaker  and 
Refinetti,  1992).  In  these  studies,  wild-type  hamsters 
were  initially  maintained  on  24-h  and  tau  mutants  on 
20-h  light-dark  (LD)  cycles.  The  animals  were  then 
transferred  to  constant  darkness  (DD)  and  PRCs  sub¬ 
sequently  measured  after  7  days  and  49  days.  When 
tested  after  7  days  in  DD,  wild-type  and  tau  mutant 
hamsters  exhibited  identical  type  1  PRCs  (Fig.  1A). 
Remarkably,  during  49  days  in  DD,  the  tau  mutant 
PRC  amplitude  increased  progressively,  until  a  transi¬ 
tion  to  type  0  PRC  was  observed  (Fig.  IB).  The 
wild-type  PRC,  on  the  other  hand,  remained  largely 
unchanged  after  49  days  in  DD  (data  not  shown). 

The  second  remarkable  behavioral  effect  of  pro¬ 
longed  exposure  to  DD  in  hamsters  was  the  gradual 
lengthening  of  the  activity  time  (a)  observed  in  both 
wild-type  (S.  Yamazaki,  data  not  shown)  and  tau 
mutant  hamsters  (Fig.  1C).  In  another  set  of  experi¬ 
ments,  groups  of  tau  mutant  hamsters  were  main¬ 
tained  in  various  T-cycles  and  then  subjected  to  light 
pulses.  In  these  experiments,  there  was  a  direct  corre¬ 
lation  between  the  amplitude  of  phase  shifts  for  each 
experimental  group  and  the  corresponding  a  lengths 
caused  by  the  different  T-cycles  (Shimomura  and 
Menaker,  1994). 

Pittendrigh  and  Daan  (1976b)  proposed  that  the 
pacemaker  in  many  animals  is  com  pound,  com  prising 
at  least  evening  (E)  and  morning  (M)  oscillators.  Two 
bouts  of  activity  are  considered,  the  first  controlled  by 
the  E  oscillator  and  the  second  by  the  M  oscillator.  The 
resulting  activity  has  duration  a;  hence,  the  value  of  a 
reflects  the  phase  relationship  between  E  and  M.  Their 
proposition  was  based  on  the  observation  of  "split¬ 
ting"  of  animal  activity  into  two  distinct  components 
and  of  bimodal  activity  peaks  in  some  animals.  These 
often-covert  components  become  evident  when  light¬ 
ing  conditions  are  altered.  This  two-oscillator  model 
was  explored  theoretically  by  Daan  and  Berde  (1978), 


Figure  1.  Phase  response  curves  measured  at  2-h  intervals.  (A) 
PRC  of  homozygous  tau  mutant  hamsters  (filled  circles,  SS)  and  a 
wild-type  PRC  (open  triangles,  WT)  obtained  by  Takahashi  et  al. 
(1984)  after  7  days  of  DD  are  replotted  on  equivalent  circadian 
coordinates.  Standard  errors  are  indicated  by  the  bars  or  are 
within  the  symbols.  (B)  PRCs  of  homozygous  tau  mutant  ham¬ 
sters  after  7  (filled  circles)  and  49  (open  circles)  days  in  constant 
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darkness.  (C)  Relationship  between  increase  in  activity  time 
(filled  circles,  alpha)  and  increase  in  light-induced  phase  delay 
(open  triangles,  phase  delay)  with  increased  number  of  days  in 
constant  darkness  in  homozygous  tau  mutant  hamsters.  Alpha 
estimated  by  naive  observers;  bars  indicate  standard  errors.  From 
Shimomura  and  Menaker  (1994).  l*AUTHORt  FL E A-S-fe 
ppnyjnr  rtrrrTrr  was  TAKFN^nrpCTLY 

FROM.  XaiS-SWRePH — 


ing  oscillators  are  identical;  thus,  subscripts  are  omit¬ 
ted  when  they  refer  to  identical  parameter  values.  The 
following  equations  describe  the  coupled  oscillator 
system: 

Evening  oscillator  (E): 

%>R£-cS£-bS£2+(d-L)  +  K,  (1) 


Kawato  and  Suzuki  (1980),  Kawato  (1985),  and  Mori  et 
al.  (1994).  The  two-oscillator  model  explains  the 
results  of  many  experiments  on  hamsters  (Pittendrigh 
and  Daan  1976b;  Elliot  and  Tamarkin,  1994;  Gorman  et 
al.,  1998),  and  therefore  we  employ  it  here. 

The  coupling  of  oscillators  adds  greatly  to  the  diffi¬ 
culty  of  understanding  the  interdependence  of  x  and 
PRC  shapes  because  it  is  necessary  to  consider  the 
intrinsic  properties  of  the  constituent  oscillators,  the 
nature  of  the  coupling,  and  the  emergent  properties  of 
the  coupled  system.  We  describe  here  our  computer 
simulations  of  the  tau  mutant  hamster  pacemaker 
viewed  as  a  system  of  mutually  coupled  E  and  M  oscil¬ 
lators.  Our  aim  was  to  determine  the  relationships 
between  oscillator  amplitude,  period,  PRC,  and  activ¬ 
ity  duration.  During  these  analyses,  we  found  that 
two  complementary  approaches  are  useful  for  analyz¬ 
ing  such  coupled  oscillator  systems.  When  the  system 
is  weakly  coupled,  we  can  employ  either  a  "distinct 
oscillators"  view  in  considering  the  properties  of  the 
separate  components  and  a  "system"  view  when  con¬ 
sidering  the  emergent  properties  of  the  system  as  a 
whole.  For  strongly  coupled  systems,  only  the  system 
view  is  appropriate. 


THE  MODEL 

Our  modeling  analyses  were  aimed  primarily 
toward  an  analysis  of  the  results  of  experiments  on 
hamsters  performed  by  Shimomura  and  Menaker 
(1994). 

Our  model  pacemaker  system  is  based  on  a  set  of 
two  identical  Pittendrigh-Pavlidis  equations  to  simu¬ 
late  the  E  and  M  oscillators.  These  equations,  devel¬ 
oped  by  Pavlidis,  capture  many  of  the  formal  proper¬ 
ties  of  the  Drosophila  eclosion  pacemaker  (Pavlidis, 
1967)  and  were  examined  extensively  by  Pittendrigh 
etal.  (1991).  Each  oscillator,  E  or  M,  is  described  by  two 
first-order  nonlinear  differential  equations.  State  vari¬ 
ables  (R  and  S,  which  are  the  same  letters  employed  by 
Pittendrigh)  are  indexed  by  E  and  M  to  designate  the 
specific  oscillator.  We  assume  that  evening  and  morn- 


‘s^  =  RE-aS£+CW£SM.  (2) 

Morning  oscillator  (M): 

dRM/«  =  RM-cSM-bS2M+(d-L)  +  K,  (3) 

dSM/dl=^M  ~a$M  W 

Four  parameters  describe  the  morning  and  evening 
oscillators:  a ,  b,  c,  and  d.  Each  set  of  equations  includes 
a  small  nonlinear  term,  K  (K  =  kj (1  +  Jc,R2),  kx  =  1,  k2  = 
100),  formulated  by  W.  T.  Kyner  (Pittendrigh  and 
Kyner,  personal  communication,  *YEAiR?*)  to  make 
the  equations  smoother.  The  effect  of  light  on  the  sys¬ 
tem  is  mediated  by  the  term  L.  Positive  values  of  L  cor¬ 
respond  to  light  pulses  of  amplitude  L.  Light-pulse 
duration  is  that  interval  during  which  L  is  nonzero. 

We  define  AM  and  AE,  respectively,  as  the  ampli¬ 
tudes  of  the  M  and  E  oscillators,  with  each  term  com¬ 
puted  as  follows  (see  appendix): 

Am  =  A(RU)  +  A(Sm)  and  A£  =  A{RL)  +  A(SC). 


i  i 
A 


The  amplitude  term,  A(RJ,  is  the  peak-to-trough 
excursion  of  state  variable  R  for  the  morning  oscillator, 
and  A(Sm)  is  the  peak-to-trough  excursion  of  the  S 
variable  for  the  morning  oscillator.  The  amplitudes  for 
the  evening  oscillators  have  similar  nomenclature 
(Fig.  2A).  a 

The  coupling  term  CME  defines  the  effect  that  tlje  M 
has  on  E;  equivalently,  CEM  defines  the  effect  that  E  has 
on  M.  Coupling  is  linear,  mediated  through  actions  on 
the  S  variable  for  each  oscillator.  For  mostsimulations, 
coupling  strengths  between  oscillators  are  of  equal 
amplitude  but  have  opposite  sign,  with  positive  CME 
and  negative  C£M.  Positive  CME  implies  that  the  morn¬ 
ing  oscillator  drives  the  evening  one.  Similarly,  nega¬ 
tive  CM£  causes  the  evening  oscillator  to  retard  the 
morning  one.  The  consequence  of  this  coupling 
scheme  is  that  E  phase-leads  and  has  smaller  ampli¬ 
tude  than  M,  even  though  the  equations  for  the  two 
oscillators  are  otherwise  identical. 
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Figure  2.  Simulations  of  the  E  and  M  oscillators.  The  graphs  are  of  CircadianDynamix-generated  plots  of  the  state  variables  (Rt,  SE,  RM,  and 
SJ  of  the  E  and  M  oscillators.  (A)  Time-series  graphs,  wild-type,  7  days  in  DD.  The  upper  and  lower  traces  show  the  simultaneous  and 
instantaneous  values  of  the  four  state  variables  of  the  pacemaker  system.  Note  that  the  E  oscillator  has  lower  amplitude,  and  phase  leads 
the  M  oscillator.  Amplitudes  (lengths  of  bars:  A(Rt),  A(St)f  A(RM),  and  AiSJ)  are  the  peak-to-trough  excursions  of  the  respective  variables. 
The  duration  of  simulated  hamster  running  wheel  activity,  a,  is  the  interval  between  the  beginning  of  activity  set  by  the  evening  oscillator 
(time  point  at  which  Rt  exceeds  2/3  of  its  maximum  value-dashed  line  in  the  upper  panel)  and  the  end  of  the  activity  set  by  the  M  oscillator 
(time  point  at  which  RM  exceeds  2/3  of  its  maximum  value-dashed  line  in  the  lower  panel).  (B)  Phase  plane  plots  of  state  variables  for 
wild-type  and  tau  mutant  hamstersimulations.  The  simulated  system  comprises  two  limit  cycle  oscillators,  each  defined  by  R  and  S  vari¬ 
ables.  Note  the  smaller  size  of  the  limit  cycles  in  the  tau  mutant  simulations.  Coupling  strengths  of  the  simulations  of  hamsters  after  7  days 
DD  were  set  to  Cut  =  0.2  and  CIM  =  -0.2.  For  the  simulations  of  hamsters  after  49  days  DD,  they  were  set  to  CMI  =  0.05  and  CEM  =  -0.05.  Parame¬ 
ter  values  for  the  wild  type  are  (a  =  0.85,  b  =  0.3,  c  =  0.8,  d  =  0.5)  and  for  the  tau  mutant  are  ( a  =  0.85,  b  =  0.3,  c  =  1.5,  d  =  0.5). 
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We  used  Pittendrigh's  standard  set  of  parameters, 
except  for  the  parameter  a,  as  follows:  a  =  0.85,  b  -  0.3,  c 
=  0.8,  d  =  0.5  (developed  for  Drosophila  eclosion  simu¬ 
lations)  (Pittendrigh  et  al.,  1991)  to  generate  a  24-h 
"wild-type"  circadian  oscillator.  The  Pittendrigh- 
Pavlidis  equations  were  used  for  their  great  conve¬ 
nience, even  though  they  lead  to  PRCs  thathave  larger 
delay  than  advance  regions,  opposite  to  what  the 
hamsters  show  (Daan  and  Pittendrigh,  1976;  Johnson, 
1992).  This  difference,  however,  affects  only  one 
aspect  of  our  work,  as  we  will  indicate  later.  To  obtain 
simulations  of  the  tau  mutant  hamster  pacemaker,  we 
altered  only  one  value  in  this  basic  parameter  set; 
namely,  we  increased  the  value  of  c  from  0.8  to  1.5  (see 
appendix).  This  single  change  reduces  the  oscillator 
period  to  20  h  and  hence  mimics  the  period  of  the 
homozygous  tau  mutant  hamster.  All  other  parame¬ 
ters  are  identical  for  both  wild-type  and  tau  mutant 
simulations.  Circadian  time  (CT)  is  normalized,  so 
that  one  circadian  hour  always  corresponds  to  the 
period  x  divided  by  24. 

In  this  model,  the  analog  of  running  wheel  activity 
occurs  every  time  the  variable  R  in  either  the  morning 
or  evening  oscillator  is  above  some  threshold  value, 
which  we  set  to  2/3  of  the  amplitude  of  this  variable 
(Fig.  2A).  As  a  result,  we  have  an  activity  band  whose 
total  duration  includes  a  firstsubband  controlled  by  E 
and  a  second  one  controlled  by  M.  The  length  of  the 
total  activity  (a)  thus  reflects  the  phase  relationship 
between  E  and  M  oscillators.  The  onset  of  activity  was 
assigned  aj  the  phase  reference  point  for  circadian 
time  ]2fffr  h.  Model  assumptions  are  summarized  in 
the  Discussion  section. 

We  explored  the  properties  of  the  coupled  pace¬ 
maker  system  described  by  these  equations  with 
CircadianDynamix,  a  computer  program  that  is  an 
extension  of  NeuroDynamix,  originally  developed  to 
explore  the  properties  of  neurons  and  small  neuronal 
networks  (Friesen  and  Friesen,  1994;  Angstadt  and 
Friesen,  1995).  (Another  variation  of  this  program, 
CalciumDynamix,  helped  provide  new  insights  into 
mechanisms  underlying  intracellular  calcium  oscilla¬ 
tions  [Friesen  et  al.,  1995].)  In  this  program,  the 
responses  of  oscillators  can  be  simulated  under  condi¬ 
tions  of  constant  darkness,  constant  light,  and  single- 
or  double-pulse  light  entrainment  paradigms.  The 
instantaneous  values  of  state  variables  are  displayed 
as  time-series  graphs  or  as  phase  plane  plots.  Simu¬ 
lated  animal  running  wheel  activity  is  shown  as 
actograms,  with  onset  of  activity  at  CT  12.  Finally, 
CircadianDynamix  displays  instantaneous  values  for 


many  computed  quantities,  including  amplitude, 
period,  and  phase.  Two  sets  of  paired  oscillators  are 
included  in  the  program.  An  "experimental"  set  of 
oscillators  can  be  manipulated  by  light  sources.  A  sec¬ 
ond  set  of  identical  oscillators,  which  are  unaffected 
by  such  manipulations,  serve  as  controls.  Thus,  the 
values  of  simulated  light-induced  phase  shifts  are 
determined  by  a  comparison  between  phase  reference 
points  on  the  experimental  and  control  oscillators. 


SIMULATIONS 

One  of  our  central  aims  was  to  provide  an  explana¬ 
tion  for  the  differences  in  the  PRCs  of  wild-type  and 
tau  mutant  hamsters  assessed  7  days  and  49  days  in 
DD  after  being  held  under  LD  14:10  and  LD  11.7:8.3 
conditions,  respectively  (Shimomura  and  Menaker, 
1994).  Our  central  assumption  is  that  prolonged  expo¬ 
sure  to  DD  reduces  the  coupling  strength  between  the 
morning  and  evening  oscillators  in  both  wild-type 
and  tau  mutant  hamsters.  The  idea  is  that  after  7  days 
in  DD,  hamsters  exhibit  the  aftereffects  of  their  previ¬ 
ous  light  schedules,  thus  showing  a  configuration  of 
stronger  coupling.  After  49  days  in  DD,  on  the  other 
hand,  they  exhibit  their  new  state,  that  is,  weaker  cou¬ 
pling.  Thus,  for  simulations  of  the  pacemakers  in  both 
hamster  strains  after  7  days  in  DD,  we  set  the  coupling 
strengths  (CM£  and  C£J  to  relatively  large  values.  To 
simulate  the  state  of  the  pacemaker  systems  after  pro¬ 
longed  DD  (49  days),  we  reduced  the  values  of  the 
coupling  constants.  The  resulting  limit  cycles  of 
wild-type  and  tau  mutant  hamsters'  dual  pacemakers 
are  shown  in  Fig.  2B.  In  all  simulations,  the  duration  of 
the  light  pulse  was  set  to  1,0  h  with  amplitude  1.1  — an 
arbitrary  unit  of  intensity. 

For  maximum  realism,  we  simulated  the  effect  of 
long-term  exposure  to  DD  by  a  stepwise  decrease  in 
the  coupling  strength  between  the  E  and  M  oscillators, 
for  both  control  and  experimental  oscillator  systems. 
With  coupling  strength  for  CM£ranging  from  0.2  to  0.05 
(correspondingly,  ranging  from  -0.2  to  -0.05),  we 
conducted  computer  experiments  to  determine  the 
shapes  and  amplitudes  of  PRCs  in  bo4h  wild-type  and 
tau  mutanthamsters.Themost  striking  results  of  these 
simulations  are  that  changes  in  coupling  strength 
hardly  altered  the  shape  of  the  PRC  for  simulations  of 
wild-type  animals  (Fig.  3A),  whereas  the  PRC  ampli¬ 
tude  in  tau  mutant  simulations  increased  progres¬ 
sively  as  coupling  strength  was  decreased  (simulating 
progressively  longer  exposure  to  DD).  At  the  weakest 
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A, 


Figure  3.  Effect  of  coupling  strength  on  PRC  amplitude.  The  graphs  illustrate  PRC  alterations  during  long-term  exposure  of  hamsters  to 
DD  simulatedasastepwise  decrease  in  the  coupling  strength  between  EandM.(A)PRCofwild-type  pacemaker  (a  =  0.85,  b  =  0.3,  c  =  0.8,  d  = 
^  0.5).  Reduction  of  coupling  strengths  has  little  if  feet  on  PRC  amplitude  or  shape.  (B)  PRC  of  tau  mutant  pacemaker  (a  =  0.85,  b  =  0.3,  c=  1.5,  d 

V  =  0.5).  Note  that  decreasing  coupling  strengths  causes  increases  in  PRC  amplitude.  For  the  weakest  coupling,  C4,  there  is  a  transition  to  a 

type  0  PRC.  Light-pulse  intensity:  L  =  1.1;  pulse  duration  =  1.0  h  (clock  time,  not  circadian  hours).  PRCs  were  obtained  with 
CircadianDynamix  programmed  to  simulate  identical  E  and  M  Pittendrigh-Pavlidis  oscillators.  Coupling  strength  [CME  =  C,  CEM  =  -Cl:  C, 
(open  circles)  =  0.25;  C,  (open  squares)  =  0.2;  C,  (open  triangles)  =  0.1;  C4  (*)  =  0.05. 


coupling  strength  employed,  the  tau  mutant  PRC 
showed  a  transition  from  type  1  to  type  0  (Fig.  3B). 
These  results  successfully  replicate  the  PRC  alter¬ 
ations  observed  in  the  experiments  on  hamsters  (Fig. 
1).  We  also  found  that  reduced  coupling  strength 
increases  the  values  of  ain  simulations  of  the  activity 
rhythms  of  both  systems. 


Two  Views  of  Coupled  Oscillators 

To  understand  the  mathematical  basis  of  the  simu¬ 
lated  behaviors  described  above,  we  first  sought  to 
understand  the  relationship  of  the  PRC  of  the  coupled 
pacemaker  system  to  those  of  the  constituent  oscilla¬ 
tors.  That  is,  we  asked,  "For  the  hamster,  is  the  system 
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PRC  some  average  of  the  intrinsic  PRCs  of  E  and  M 
oscillators?"  If  so,  how  is  the  system  PRC  amplitude 
and  shape  affected  by  the  phase  relationship  between 
the  unit  oscillators?  Here  we  argue  that  the  coupled 
system  should  be  viewed  differently,  depending  on 
the  strength  of  coupling  between  the  oscillator  units. 
Most  generally,  the  coupled  system  has  to  be  consid¬ 
ered  as  a  whole  unit,  and  we  cannot  associate  the  PRC 
of  the  system  with  the  intrinsic  PRCs  of  the  compo¬ 
nents.  We  call  this  a  "system"  view.  Nevertheless, 
weak  coupling  strength  ensures  preservation  of  the 
intrinsic  properties  for  each  oscillator.  In  this  case,  we 
can  estimate  the  PRC  of  the  system  from  the  intrinsic 
PRC  of  the  components.  For  such  weak  coupling,  we 
can  adopta  distinct  oscillator  view  for  some  analytical 
purposes.  We  consider  the  two  views  in  this  section. 
The  experiments  by  Shimomura  and  Menaker  (1994) 
can  only  be  explained  by  adopting  the  system  view, 
but  some  features  of  the  hamsters'  oscillators  are 
better  understood  by  adopting  the  distinct  oscillators 
view. 

The  System  View 

When  the  coupling  is  strong,  the  unit  E  and  M  oscil¬ 
lators  begin  to  lose  their  original  identities,  including 
intrinsic  periods,  amplitudes,  wave  shapes,  and  PRCs. 
Instead,  the  properties  meld  into  the  corresponding 
properties  of  the  system.  Therefore,  we  cannot  esti¬ 
mate  the  PRC  of  the  coupled  system  from  the  intrinsic 
PRCs  of  free  oscillators.  Instead,  the  PRC  can  only  be 
described  for  the  system.  Furthermore,  with  strong 
coupling,  the  activity  bands  controlled  by  E  and  M 
become  cohesive  because  the  phase  relationship 
between  them  is  always  maintained.  That  is,  the 
beginning  and  end  of  the  activity  bands  tend  to  shift 
equally  to  the  steady-state  value  when  subjected  to 
light  pulses. 

According  to  Lakin-Thomas  et  al.  (1991),  the  phase 
shift  magnitude  of  an  independent  oscillator  depends 
on  (1)  impulse  intensity  (the  greater  the  impulse,  the 
larger  the  phase  shift)  and  (2)  the  amplitude  of  oscilla¬ 
tions  (i.e.,  smaller  oscillator  amplitudes  lead  to  the 
larger  phase  shifts).  Here  we  extend  Lakin-Thomas's 
principle  to  the  coupled  oscillator  system.  We  first 
define  the  amplitude  for  a  coupled  system  comprising 
two  unit  oscillators  and  then  analyze  the  interrelation¬ 
ships  between  systems  properties,  including  coupling 
strength,  a,  PRC  amplitude,  and  system  amplitude. 

We  completed  a  series  of  three  simulations  to  study 
the  nature  of  strongly  coupled  systems  using 


Table  1.  Summary  of  results  from  simulation  set  1.  Arrows  indi¬ 
cate  changes  in  parameter  values  (column  1)  and  the  consequences 
for  the  system  (columns  2-4). 
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CircadianDynamix.  Our  principal  aim  was  to  under¬ 
stand,  systematically,  why  a  direct  correlation  exists 
between  a  (which  reflects  the  phase  relationship 
between  E  and  M)  and  phase  shift  magnitude  in  tau 
mutant  hamsters  (Fig.  1C).  We  show  below  that  this 
correlation  is  indeed  mediated  by  the  system  ampli¬ 
tude  changes  as  follows.  Through  simulation  set  1,  we 
show  that  phase  shift  magnitude  is  determined 
directly  by  the  component  amplitudes  of  the  system, 
Ae  and  Am,  and  not  by  a.  In  simulation  set  2,  we  show 
that  this  magnitude  is  indeed  determined  by  the  sum 
of  A£and  AM,  which  we  call  total  amplitude  (Ar)  of  the 
system .  Finally,  via  simulation  set  3,  we  show  that  the 
negative  correlation  between  PRC  and  ATholds  gener¬ 
ally;  it  is  independent  of  specific  system  parameters. 

Simulation  Set  1:  Interrelationships  between  PRC 
Amplitude,  a,  and  Oscillator  Amplitudes 

In  these  simulations,  we  determined  if  phase  shift 
magnitude  is  affected  directly  by  a  or  indirectly,  via 
oscillator  amplitudes  A£  and  AM.  Our  approach  was 
systematically  to  vary  oscillator  parameters  while 
maintaining  a  constant  coupling  strength.  The  four 
system  parameters — a,  b,  c,  and  d — that  define  E  and 
M  oscillator  properties  provide  a  rich  repertoire  of 
possible  simulations.  Varying  these  parameters  in 
both  E  and  M  can  generate  positively  and  negatively 
correlated  changes  in  x,  a,  AE,  and  AM. 

We  first  set  up  a  basic  system  configuration  defined 
by  the  nominal  parameter  set  a  -  0.85,  b  =  0.3,  c  =  0.8, 
and  d  =  0.5  with  moderately  weak  coupling  strength 
CM£=  -C^  =  0.05.  After  recording  the  values  of  x,  a,  A M, 
and  A £,  we  changed  one  parameter  at  a  time,  while 
maintaining  all  the  others  fixed,  and  noted  the  new 
values  of  these  four  quantities.  A  sample  result  is  that 
selectively  increasing  parameter  a  leads  to  an  increase 
in  t,  increases  in  both  A £  and  A M,  and  a  decrease  in  a.  In 
Table  1,  we  present  a  qualitative  overview  of  our 
results,  demonstrating  the  interrelationships  between 
changes  in  system  parameters  and  x,  A£and  AM,and  a. 

Table  1  shows  that  amplitudes  A£and  Aw  like  a,  are 
negatively  correlated  with  parameter  b.  For  changes  in 
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Figure  4.  PRCs  of  coupled  oscillators  when  system  parameters  are  altered, changing  amplitudes  and  phase  relationships  between  E  and  M 
oscillators.  (A)  Series  of  PRCs  generated  by  decreasing  parameter  d.  These  sequential  parameter  changes  (1,  2,  and  3)  led  to  increases  in 
PRC  amplitude  and  simultaneous  decreases  in  and  AM,  while  a  increased  (Al  >  A2>  A3;  al  <  a 2  <  a3).  (B)  Series  of  PRCs  generated  by 
increasing  parameter  b.  These  parameter  changes  (1,2,  and  3)  led  to  increases  in  PRC  amplitude  and  decreases  in  AE,AM,  and  in  a(Al  >  A2> 
A3;  al  >  al  >  oc3).  Parametervalues  in  A  ( a  =  0.85;  b  =  0.3;c  =  0.8;  d  =  0.5, 0.3,  and  0.1)  and  in  B  ( a  =  0.85;  c  ~  0.8;  d  =  0.5;  b  —  0.3,  0.5,  and  1.0).  Cou¬ 
pling  strength  was  constant  for  all  simulations:  CMl  =  0.05;  CtM  =  -0.05;  pulse  intensity  =  1.1;  pulse  duration  =  1  h. 


parameter  d,  the  quantities  have  opposite  correlations; 
namely,  the  amplitudes  are  positively  correlated  with 
d,  but  for  the  activity,  the  correlation  is  negative.  How 
does  the  PRC  amplitude  correlate  with  b  and  d ?  We 
obtained  quantitative  measures  of  PRCs  for  three 
decreasing  values  of  parameter  d  (leading  to  increas¬ 
ing  a  and  decreasing  AE  and  AJ.  As  shown  in  Fig.4A, 
these  stepwise  decreases  in  d  lead  to  a  progressive 
increase  in  PRC  amplitude.  Similar  analysis  of  the 
effects  of  increasing  values  of  b  revealed  that  PRC 
amplitude  increased  progressively  even  with  a 
decreasing  aand  decreasing  A£and  AM(Fig.  4B).  These 
simulations  demonstrate  that  the  PRC  amplitude  is 
negatively  correlated  with  oscillator  amplitude  but  is 
not  necessarily  positively  or  negatively  correlated 
with  values  for  a.  We  generalize  this  result  below. 

Simulation  Set  2:  Interrelationships  between  Cou¬ 
pling  Strength,  a,  and  System  Amplitude 

In  these  simulations,  we  define  the  total  amplitude 
A,  of  the  coupled  system  and  show  that  PRC  ampli¬ 
tude  is  negatively  correlated  to  Ar  Having  described 
the  interrelationships  of  amplitude,  activity  duration, 
and  PRC  amplitude  when  AE  and  AM  covary,  we  now 
explore  systems  showing  increasing  AE  and  decreas¬ 
ing  Am.  Such  systems  are  obtained  by  setting  oscillator 


parameters  to  a  fixed  value  and  then  altering  only  the 
coupling  strength.  As  in  the  simulations  described 
above,  CME and  C£M,  parameters  that  describe  the  effect 
ofM  on  E  and  of  E  on  M,  respectively,  always  have  the 
same  amplitudes  but  opposite  sign,  with  CM£  positive 
and  C£M  negative. 

We  found  that  decreasing  coupling  strength  always 
leads  to  increases  in  a,  that  is,  to  a  lengthening  of  the 
simulated  activity.  Stated  in  other  words,  decreasing 
the  coupling  strength  increases  the  phase  difference 
between  E  and  M.  Another  consequence  of  reducing 
coupling  strength  is  that  AE  increases  and  AM 
decreases,  a  direct  consequence  of  positive  CM£  and 
negative  Cm.  To  carry  out  more  quantitative  analyses, 
we  defined  the  total  amplitude,  Av  of  the  coupled  sys¬ 
tem  as  the  sum  of  the  E  and  M  oscillator  variable 
amplitudes: 

Ar  =  A£  +  AM  -  A(RC)  +  A(Sf)  +  A(RJ  +  A(SJ. 

We  determined  the  explicit  values  assumed  by  AE,  AM, 
and  AT  for  a  wide  range  of  coupling  strengths,  span¬ 
ning  those  we  employed  to  simulate  brief  and  pro¬ 
longed  exposure  to  DD  in  wild-type  (Fig.  5A)  and  tau 
mutant  hamsters  (Fig.  5B).  The  PRC  amplitudes 
depicted  in  Fig.  3  can  now  be  understood  in  light  of  the 
relationship  between  oscillator  amplitude  and  cou- 
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Figure  5.  Comparison  of  am  plitude  changes  in  E  and  M  oscillators  with  those  of  the  coupled  system  as  coupling  strength  is  decreased.  (A) 
Simulations  of  the  wild-type  pacemaker.  (B)  Simulations  of  the  tau  mutant  pacemaker.  The  abscissa  in  each  graph  is  the  size  of  the  cou¬ 
pling  coefficient— positive  for  CMland  negative  for  CEM.The  insets  describe  the  values  for  coupling  strength  employed  for  the  simulations 
illustrated  in  Figure  3.  Oscillations  dam  p  out,  no  longer  self-sustained,  for  coupling  strengths  greater  than  about  0.39.  System  parameters 
were  those  described  in  the  caption  of  Figure  3. 

pling  strength.  In  simulations  of  wild-type  hamsters  found  that  as  coupling  strength  decreases,  the  ampli- 

(Fig.  5A),  we  found  that  as  coupling  strength  tude  of  E  increases,  the  amplitude  of  M  decreases,  and 

decreases,  the  amplitude  of  E  (A £)  increases,  the  am pli-  the  total  amplitude  also  decreases.  These  substantial 

tude  of  M  (Am)  decreases,  and  the  total  amplitude  AT  changes  in  ATare  reflected  in  the  large  changes  in  PRC 

remains  nearly  constant.  Please  note  again  that  the  amplitude,  which  are  negatively  correlated  to  AT( cf. 

PRCs  generated  for  these  different  coupling  strengths  Figs.  3B  and  5B). 

(Fig.  3A)  are  all  of  nearly  equal  amplitude.  On  the  For  the  wild-type  system,  AE  and  AM  changes  have 

other  hand,  for  the  tau  mutant  system  (Fig.  5B),  we  equal  absolute  values  but  are  of  opposite  signs.  Thus, 
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their  sum,  Av  is  nearly  unaltered,  and  consequently 
PRC  amplitudes  remain  nearly  unchanged.  On  the 
other  hand,  in  the  tau  mutant  system,  AM  decreases 
much  more  than  AE  increases;  thus,  A T  decreases  with 
exposure  to  DD.  Therefore,  the  tau  mutant  PRC 
increases  continuously  during  the  decreasing  cou¬ 
pling  process  associated  with  prolonged  exposure  to 
DD.  With  sufficiently  long  exposure  to  DD,  the  tau 
mutant  PRC  changes  to  type  0.  The  reasons  for  this 
transition  are  described  below  in  the  final  simulation 
set. 

Simulation  Set  3:  Transitions  between 
Type  1  and  Type  0  PRCs 

Correlations  between  oscillator  properties,  such  as 
the  value  of  aand  system  amplitude,  usually  depend 
on  the  choice  of  the  parameter  sets.  To  demonstrate  the 
generality  of  the  behaviors  observed  in  the  simulation 
sets  described  above,  we  tested  the  following  proposi¬ 
tion:  if  the  phase  shift  magnitude  depends  only  on  Av 
then  there  must  be  a  threshold  value  of  Arat  which  a 
transition  from  type  1  to  type  0  PRC  occurs,  independ¬ 
ent  of  a  and  of  the  parameter  sets  used. 

We  determined  PRCs  for  coupled  oscillators 
defined  by  a  broad  range  of  parameter  values  ( afb,c,d , 
CtM,  and  CAJ.  Each  set  of  parameters  defines  a  differ¬ 
ent  configuration  for  the  coupled  oscillator  system 
and  generates  a  unique  pair  of  values  for  A  T  and  a  (see 
Fig.  6).  In  this  illustration,  parameter  sets  that  generate 
type  1  PRCs  are  represented  by  open  symbols,  and 
those  that  generate  type  0  PRCs  are  denoted  by  filled 
symbols.  The  plot  shows  clearly  that  the  zone  in  which 
transitions  from  type  1  (open  symbols)  to  type  0  (filled 
symbols)  occurs  is  independent  of  a.  Such  transitions 
are  confined,  moreover,  to  a  very  narrow  band  of  val¬ 
ues  for  Ar  Thus,  systems  characterized  by  large  values 
of  A,  have  type  1  PRCs,  and  those  with  low  total 
amplitudes  have  type  0  responses.  We  also  demon¬ 
strated  that  this  is  true  for  systems  presenting  differ¬ 
ent  coupling  signs,  such  as  (+CME,  +CEM)  and  (-CME, 
-CJ  (data  not  shown).  System  configurations  dis¬ 
playing  equaM7but  very  different  a  values  have  iden¬ 
tical  PRCs  (data  not  shown).  As  is  well  known,  PRC 
amplitude  depends  on  the  intensity  of  the  applied 
light  stimulus  light  (Winfree  1980);  therefore,  the  zone 
for  type  1  to  type  0  transitions  is  different  for  every 
light  impulse  intensity.  In  particular,  as  the  intensity  of 
the  light  pulse  is  raised,  these  transitions  occur  at 
greater  AT  values  (data  not  shown). 


Remembering  that  system  configurations  are  set 
both  by  parameters  of  the  unit  oscillators  and  by  cou¬ 
pling  coefficients,  we  link  systems  that  have  identical 
component  oscillators  but  with  decreasing  coupling 
strengths  (arrows,  Fig.  6).  We  also  indicate  the  average 
range  of  Tin  these  linked  systems.  These  data  demon¬ 
strate  that  the  value  of  a  consistently  increases  as  cou¬ 
pling  strength  decreases  (direction  of  the  arrows)  and 
that  At  can  increase  or  decrease  with  changes  in  cou¬ 
pling  strength.  Thus,  for  the  wild-type  pacemaker  (f  = 
24  h),  decreasing  the  coupling  strength  causes  only  a 
small  decrease  in  Av  and  these  are  far  from  the  transi¬ 
tion  zone.  On  the  other  hand,  decreasing  coupling 
strength  causes  larger  decreases  in  Ar  of  the  tau 
mutant  hamster  system  (f  =  20  h),  and  these  values  of 
ATare  near  the  transition  zone.  Hence,  the  tau  mutant 
system  makes  the  transition  to  a  type  0  PRC  because 
its  amplitude  is  more  sensitive  to  reductions  in  cou¬ 
pling  strength  and  because  this  pacemaker  system  is 
adjacent  to  the  type  1  to  type  0  transition  zone. 

Distinct  Oscillator  View 

The  system  view  must  be  adopted  at  any  coupling 
strength  to  analyze  interrelationships  between  x,  a, 
and  amplitudes.  As  coupling  strength  is  weakened, 
however,  the  influence  of  each  component  oscillator 
on  the  dynamics  of  the  other  is  weakened,  allowing 
each  component  to  gradually  recover  its  intrinsic 
properties,  even  if  not  completely.  In  such  a  weakly 
coupled  system,  the  properties  of  the  individual  oscil¬ 
lators  become  evident  and  accessible  to  analysis.  The 
distinct  oscillators  view  is  well  suited  to  explain  com¬ 
mon  phenomena  of  hamster  circadian  pacemaker, 
such  as  the  differential  shifts  of  the  beginning  and  end 
of  the  activity  bands  in  response  to  light  pulses. 

Light-induced  phase  shifts  of  the  weakly  coupled  E 
and  M  pacemaker  system  reflect  the  combined,  intrin¬ 
sic  phase  response  curves  of  E  and  M  oscillators  (Mori 
et  al.,  1994).  Two  PRCs  for  a  definite  light-pulse  inten¬ 
sity,  PRC£  and  PRCM,  generated  from  simulations  of 
uncoupled  wild-type  E  and  M  oscillators,  are  depicted 
in  Fig.  7A.  We  then  applied  extremely  weak  coupling 
(C  =  0.01)  to  this  system  and  generated  a  third 
PRC— the  system  PRCEJVJ  (Fig.  7A).  In  our  simulations, 
we  are  using  two  identical  E  and  M  oscillators,  which 
are  sufficient  to  describe  the  following  behaviors.  In 
this  weakly  coupled  system,  M  lagged  £  by  3.8  h. 
Hence,  in  superimposing  the  three  PRCs  in  Fig.  7A,  we 
shifted  the  abscissa  for  PRCM  by  +3.8  h.  The  system 
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Figure  6.  Total  amplitude  (AT)  and  activity  length  (a)  of  coupled  oscillator  configurations  generated  by  a  wide  range  of  parameter  values. 
Each  point  was  generated  by  a  different  set  of  values  for  parameters  a,  b,c,  and  d  and  of  the  coupling  coefficients  CMIand  CEM  and  placed  on 
the  graph  at  the  system  amplitude  (AT)  and  activity  duration  (a).  Depending  on  the  shape  of  the  system  PRC,  we  represent  the  configuration 
as  filled  circles  if  the  PRC  is  type  1  or  as  open  circles  if  the  PRC  is  type  0.  The  symbol  A  denotes  systems  that  lie  exactly  on  this  transition, 
that  is,  when  the  pulse  phase-shifts  the  oscillator  to  a  singularity  point.  Note  that  all  type  1  systems  are  in  the  upper  region  with  high  AT  and 
that  the  type  0  region  is  restricted  to  small  values  of  A  T,  independent  of  the  a  values.  Arrows  link  systems  with  identical  values  of  a,b,c,  and 
d  but  with  coupling  strength  decreasing  in  the  direction  of  the  arrows.  As  coupling  is  weakened,  a  always  increases,  but  AT  can  either 
increase  or  decrease.  Average  I  values  are  indicated  for  four  series  of  simulations.  Note  that  the  tau  mutant  series  is  at  the  PRC  transition 
zone,  whereas  the  wild-type  series  is  far  above  this  zone. 


PRCMf  then  lies  between  those  of  E  and  M.  We  must 
emphasize  that  this  coupled  system  PRC  estimation 
from  the  intrinsic  PRCs  of  the  components  is  only 
valid  for  very  weak  coupling. 

By  convention,  the  onset  of  running  wheel  activity 
of  hamsters  is  assigned  the  phase  reference  point  CT 
12.  In  accordance  with  this  convention,  we  assigned 
the  phase  CT  12  to  activity  onset  in  the  evening  oscilla¬ 
tor  in  our  simulations  of  both  coupled  and  independ¬ 
ent  oscillators.  Thus,  CT  12  is  the  identical  reference 
point  for  the  evening  oscillator  and  for  the  coupled 
system.  For  the  specific  modeling  experiment 
described  here,a  response  to  a  lightpulse  presented  at 
phase  0  for  the  coupled  system  is  generated  by  the 
combined  responses  of  the  evening  oscillator  at  phase 
and  of  the  morning  oscillator  at  phase  <t>-  3.8  h.  As 
illustrated  in  Fig.  7 A,  there  are  two  system  phases, 


and  02,  where  there  is  an  intersection  of  PRCe  and 
PRCm.  Light  pulses  presented  at  these  phases  cause 
identical  shifts  in  E  and  M  and  hence  also  in  the  cou¬ 
pled  system.  There  are  two  phase  ranges,  ^ and  r2,  dur¬ 
ing  which  PRC£  and  PRCM  have  the  opposite  sign;  at 
all  other  phases,  PRC£  and  PRCM  have  same  sign  but 
different  amplitudes.  As  a  consequence  of  weak  cou¬ 
pling,  there  are  two  primary  properties  of  this  pace¬ 
maker  system. 

1.  Immediate  responses  of  the  E  and  M  oscillators  to  a 
pulse  of  light:  the  oscillators  undergo  nearly  inde¬ 
pendent  phase  shifts  (in  direction  and  amplitude). 
These  initial  phase  shifts  correspond  to  those 
described  by  PRC£and  PRCM  (Fig.  7A).  If  the  pulse  is 
given  at  <J>i  or  O2,  both  oscillators  undergo  the  same 
phase  shift  with  no  transient  change  in  a.  If  the  pulse 
falls  on  ranges  r\  or  X2,  one  oscillator  will  advance, 
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Figure  7.  Phase  response  curves  in  weakly  coupled  pacemaker  systems:  wild  type.  (A)  Computer-generated  PRCs  for  uncoupled  inde¬ 
pendent  oscillators,  PRCt  (filled  circles)  and  PRCM  (open  circles), and  for  the  very  weakly  coupled  system  (tosimulate  prolonged  exposure 
to  DD),  PRCtM  C).  The  abscissa  forPRCM  is  shifted  by  3.8  h. The  intrinsic  PRCs  of  £  and  M  intersect  at  <t>,  and  <t>y  Phase  ranges  r,  and  r3  denote 
phases  for  which  PRCt  and  PRCM  have  opposite  signs.  At  all  other  phases,  these  curves  have  the  same  sign  but  differing  amplitudes.  (B) 
Raster  plots  showing  simulated  actograms  of  hamsters.  Light  pulses  are  set  by  the  system  on  the  days  indicated  by  (ordinate)  atthe  fol¬ 
lowing  phases:  (I)  CT  8  (indicated  by  an  arrow  in  Fig.  7A),  (II)  CT  11.25  (at  <!>,), and  (III)  CT 17  (indicated  by  an  arrow  in  Fig.  7A, during  phase 
range  r3).  Because  of  artifacts  generated  by  the  com  pu ter  program  on  the  "day"  that  the  light  pulse  is  applied,  the  activity  for  that  day  is  not 
shown.  AtCT8  (CT  relative  to  the  £  oscillator),  both  £  and  M  are  immediately  delayed,  but  £  displays  a  much  larger  delay  than  M,  causing  a 
transient  shortening  of  a  The  £  oscillator  then  slowly  advances.  At  O,,  both  oscillators  phase-shift  equally  (II).  At  a  phase  in  r„  £  initially 
advances,  while  M  initially  delays  causing  a  transient  increase  in  a  (III).  During  subsequent  cycles,  a  returns  to  its  previous  value  and  a 
new,  steady-state  phase  is  established.  Each  line  represents  1  day,  and  the  full-scale  abscissa  is  24  circadian  hours.  (C)  Actogram  forham- 
sters  (Shimomura,  1998),  showing  a  transient  compression  in  a  due  to  a  pulse  given  at  CT  16,  in  wild-type  hamster.  Simulation  results  were 
for  the  wild -type  system  ( a  -  0.85,  b  =  0.3, c  =  0.8, d  =  0.5)  with  very  weak  coupling  (CME  =  0.01,  ClM  =  -0.01).  Data  were  generated  by  simulating 
1-h  light  pulses  with  intensity  L  =  1.1  and  duration  1.0  h.  The  PRCs  of  the  intrinsic  £  and  M  oscillators  and  for  the  coupled  system  are  type  1. 


while  the  other  will  be  delayed.  Finally,  if  the  light 
pulse  occurs  at  any  other  phase,  both  oscillators  will 
phase-shift  in  the  same  direction  but  with  differing 
initial  amplitude. 

2.  Long-term  responses  of  the  oscillators:  following  the 
initial,  transient  responses,  the  phase  relationship 
between  £  and  M  returns  to  its  original  value,  and  the 


system  phase  approaches  a  new  steady  value.  This 
system  shift  in  phase  is  described  by  PRCem  (Fig.  7A). 
The  time  interval  required  for  the  return  to  steady 
state  gives  rise  to  transients  in  phase  and  a. 
Depending  on  the  phase  of  the  light  pulse,  a  may  be 
transiently  decreased  (Fig.  7B-I)  or  increased  (Fig. 
7B-III)  for  several  cycles. 
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Such  independent  shifts  of  evening  and  morning 
oscillators  in  response  to  light  pulses  were  observed 
during  experiments  by  Elliot  and  Tamarkin  (1994)  and 
Shimomura  (1998)  (Fig.  7C).  Similarly,  we  found  that  a 
light  pulse  at  CT  8  generated  a  large  initial  delay  in  E 
and  a  smaller  delay  in  M  (Fig.  7B-I).  Subsequently,  E 
transiently  advances  and  M  delays,  reestablishing 
their  initial  phase  relationship.  This  pattern  corre¬ 
sponds  to  the  experimental  results  obtained  for  ham¬ 
sters  subjected  to  light  pulses  at  CT  16  (Fig.  7C).  For 
the  animal  experiments,  however,  the  behavior  was 
observed  in  the  advance  region,  rather  than  in  the 
delay  region,  of  the  PRC.  To  simulate  this  behavior  at 
CT  16  clearly,  we  need  large  advance  shifts  for  E  and 
M.  We  were  unable  to  simulate  this  because  our  model 
oscillator  (devised  for  Drosophila)  displays  a  PRC  with 
larger  delay  and  smaller  advance  regions  when  com¬ 
pared  to  the  PRCs  of  hamsters.  The  underlying  mech- 
anisms  that  generate  the  transient  advances  and 
delays  are  nevertheless  identical. 

The  same  relative  dynamical  independence  of  E 
and  M  oscillators  under  weak  coupling  is  even  more 
evident  for  the  tau  mutant  simulations  in  which  both  E 
and  M  have  type  0  PRCs  under  standard  stimulation 
(Fig.  8A).  Note  that  PRC£M  again  lies  between  PRCf 
and  PRCM  but  not  at  the  midpoint.  In  this  case,  PRC£ 
and  PRCM  do  not  intersect,  but  there  is  a  phase  range, 
p  r),  where  they  have  opposite  signs.  For  this  phase 

3  range,  the  activity  band  controlled  by  E  phase-delays, 

while  that  controlled  by  M  phase-advances;  thus, 
th  ese  bands  cross  (Fig.  8B).  Concomitant  with  these 
phase  shifts,  there  is  a  large,  transient  reduction  in  ocas 
the  activity  bands  cross.  A  simulation  of  the  same  sys¬ 
tem  with  strong  coupling  is  shown  in  Fig.  8C  for  com¬ 
parison.  In  this  case,  the  activity  bands  controlled  by  E 
and  M  shift  in  the  same  direction,  nearly  preserving 
the  phase  relationship  between  them. 

The  transients  in  our  simulations  are  conceptually 
different  from  those  resulting  from  a  master-slave 
oscillator  structure  (Pittendrigh,  1967).  Double-pulse 
experiments  in  Drosophila  showed  that  the  clock 
phase-shifts  immediately,  whereas  the  overt  rhythm 
shifts  slowly.  Transients  in  our  simulations  corre- 
'f^u\  K/  sPond  t0  the  time  interval  required  by  tWe  components 
j  within  the  master  clock  to  reestablish  meir  phase  reta¬ 
ke  0  tionship  after  a  light  pulse.  Elliot  and  Pittendrigh 

(abstract  from  the  1996  SRBR  meeting)  have  shown 
that  in  double-pulse  experiments  with  hamsters,  the 
clock  takes  several  cycles  to  phase-shift.  According  to 
our  model,  although  each  component  E  and  M 


phase-shifts  immediately,  the  whole  system  takes 
many  cycles  to  attain  its  new  phase. 

Large  transients  are  the  hallmark  of  weakly  cou¬ 
pled  oscillators.  Accordingly,  double-pulse  experi¬ 
ments  on  hamsters  maintained  for  different  durations 
in  DD  can  test  the  hypothesis  of  weakening  the  cou¬ 
pling  between  E  and  M  under  DD.  Under  this  hypoth¬ 
esis,  the  longer  the  animal  is  in  DD,  the  longer  the  tran¬ 
sients  within  the  clock  should  be. 

The  distinct  oscillators  view  explains  the  transient, 
differential  phase  shifts  of  the  onset  and  end  of  the 
activity  bands  in  hamsters.  These  phenomena  have 
strongly  supported  the  two-oscillator  model,  in  which 
activity  onset  and  end  are  assumed  to  be  controlled, 
respectively,  by  E  and  M  oscillators. 


DISCUSSION 


Summary  of  Results 

Physiological  experiments  provide  insights  into 
the  properties  of  specific  systems  and  their  responses 
to  perturbations.  Modeling  studies  are  valuable 
because  they  can  provide  insights  into  the  mechanistic 
origins  of  experimental  results.  More  than  that,  the 
many  configurations  generated  during  modeling 
analyses  can  give  insight  into  the  general  underlying 
principles  of  systems  thatare  independent  of  any  par¬ 
ticular  configuration.  We  have  presented  the  results  of 
experiments  using  computer  simulations  to  analyze 
the  properties  that  emerge  from  coupling  between 
two  oscillators.  Our  primary  aim  was  to  understand 
the  bases  for  alterations  in  the  PRC  amplitudes 
observed  in  wild-type  and  tau  mutant  hamsters 
exposed  to  long-term  DD.  Our  simulation  experi¬ 
ments  demonstrate  the  following: 

1.  A  change  in  a  single-system  parameter  can  mimic  the 
reduced  period  of  the  tau  mutation,  but  a  necessary 
change  occurs  also  in  other  system  properties.  In  our 
simulation,  the  parameter  change  also  decreases  the 
total  amplitude  of  the  coupled  oscillator  system. 

2.  Incorporating  the  assumption  that  the  primary  effect 
of  exposure  to  DD  is  the  gradual  weakening  of  the 
coupling  between  the  E  and  M  oscillators  that  com¬ 
prise  the  pacemaker  led  to  a  large  increase  in  PRC 
amplitude  in  oursimulated  tau  mutantsystem  butnot 
the  wild-type  system. 

3.  Although  a  direct  correlation  was  shown  between  a 
and  phase  shift  magnitude,  there  is  no  causal  relation¬ 
ship  between  them.  In  our  simulations,  as  coupling 
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Figure  8.  Phase  response  curves  in  coupled  pacemaker  systems:  tau  mutation.  (A)  Computer-generated  PRCs  for  the  uncoupled  unit  oscil¬ 
lators,  PRC,  (filled  circles)  and  PRCM  (open  circles),  and  for  the  weakly  coupled  system  (to  simulate  prolonged  DD  exposure),  PRCEM  (*). 
Note  that  all  three  PRCs  are  type  0.  In  the  coupled  system,  M  lags  £  by  5.25  h.  The  PRCs  for  £  and  M  never  cross,  but  there  is  a  range  of 
phases,  rt,  during  which  PRCt  and  PRCM  have  opposite  signs.  (B)  System  response  to  pulse  present  during  range  r,  (CT  17).  Note  that  this 
stimulus  shifts  £  and  M  activity  bands  strongly  in  opposite  directions  so  that  these  bands  cross.  (C)  Phase  shifts  in  strongly  coupled  sys¬ 
tems.  Although  the  simulated  light  pulse  was  presented  at  CT  16.5,  during  strong  coupling  causes  £  and  M  activity  bands  to  shift 
togetheras  a  unit.  (D)  Comparison  with  experimental  data  obtained  in  a  tau  mutant  hamster  (from  Menakerand  Refinetti,  1992),  showing 
the  pattern  simulated  in  B.  Oscillators  parameters  for  all  simulations  were  set  to  simulate  tau  mutant  values:  a  =  0.85,  b  =  0.3,  c  =  1.5,  and  d  = 
0.5.  The  weak  coupling  in  A  and  B  were  simulated  with  CUI  =  0.01  and  C£M  =  -0.01,  with  light-pulse  amplitude  =  1.1  and  duration  =  1  h.In  part 
C,  CM(  =  0.38,  and  CtM  =  -0.38,  with  light-pulse  amplitude  =  1.3  and  duration  =  1  h. 


strength  decreases,  a  increases  due  to  a  greater  phase 
difference  between  £  and  M.  At  the  same  time,  A j 
decreases,  causing  an  increase  in  observable  PRC 
amplitudes. 

4.  The  high  sensitivity  of  the  tau  mutant  hamster  and 
low  sensitivity  of  the  wild-type  hamster  to  prolonged 


exposure  to  DD  result  from  the  contrasting  location  of 
these  systems  on  the  total  amplitude  scale,  shown  in 
Fig.  6;  the  tau  mutant  system  total  amplitude  lies  near 
the  type  1-type  0  transition  value,  whereas  the 
wild-type  system  lies  far  from  this  value. 
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Model  Assumptions 

Several  specific  assumptions  are  incorporated  into 
our  model.  Here  we  discuss  why  these  assumptions 
are  necessary  and  how  they  affect  the  results  of  this 
analysis. 

A  distinction  should  be  made  between  the  model¬ 
ing  style  we  use  here— namely,  the  qualitative  analy¬ 
sis  of  formal  limit  cycle  oscillators  and  the  use  of  dif¬ 
ferential  equations  to  represent  explicitly  known 
cellular  processes  (Goldbeter,  1996).  In  the  first 
approach,  the  state  variables  (R  and  S  in  the  case  of  the 
Pittendrigh-Pavlidis  equation)  and  parameters  ( a ,  b,  c, 
and  d)  have  no  physiological  meaning.  They  are  sim¬ 
ply  the  components  of  the  formal  system  of  equations. 
The  aim  of  this  formal  approach  is  to  explain  qualita¬ 
tive  features  of  rhythms  in  terms  of  generic  limit  cycle 
oscillator  properties.  In  the  second  approach,  state 
variables  and  parameters  correspond  explicitly  to 
concentrations  of  chemical  substances  and  kinetic 
constants.  They  are  meant  to  provide  a  quantitative 
view  of  the  limit  cycle  oscillation  generation  and 
dynamic  properties  of  specific,  identified  circadian 
pacemakers  in  cells,  plants,  or  animals.  We  employed 
the  first  approach  here  because  the  phenomena  we 
wanted  to  investigate  can  be  explained  in  terms  of 
generic  limit  cycle  oscillator  properties.  In  addition, 
identities  of  the  state  variables  and  of  system  parame¬ 
ters  are  unknown  in  the  hamster  circadian  pacemaker. 
Either  approach  seeks  to  model  the  system  with  the 
smallest  set  of  assumptions  to  isolate  the  critical  com¬ 
ponents  and  to  focus  on  the  main  features  of  the  phe¬ 
nomena  under  investigation.  Further  refinement  of 
the  model,  with  additional  assumptions, can  be  imple¬ 
mented  if  the  initial,  simple  set  of  assumptions  cannot 
lead  to  an  explanation  of  the  phenomena.  All  the 
assumptions  shown  below  were  thus  made  primarily 
for  the  sake  of  simplicity. 

Evening  and  morning  oscillators  can  be  represented  by 
two-dimensional  limit  cycles .  Biochemical  processes 
generating  the  circadian  oscillation  probably  com¬ 
prise  many  more  than  two  proteins  and  mRNAs 
whose  concentrations  are  the  likely  state  variables  of 
cellular  circadian  oscillators  (Leloup  and  Goldbeter, 
1998;  Sangoram  et  al.,  1998).  On  the  other  hand,  we  are 
deducing  the  formal  properties  of  the  pacemaker 
through  the  direct  observation  of  the  output  rhythm, 
which  reflects  the  emergent  oscillator  properties  of  the 
complex  biochemical  processes.  This  formal  approach 
allows  us  to  describe  the  pacemaker  by  a  two-dimen¬ 
sional  minimal  representation  that  provides  phase 


and  amplitude  information  of  the  limit  cycle  and 
accounts  for  both  type  1  and  type  0  resetting  (Winfree, 
1980).  To  simulate  the  activity  controlled  by  this  pace¬ 
maker,  we  define  "activity"  to  occur  when  the  R  vari¬ 
able  of  either  evening  or  morning  oscillator  is  higher 
than  a  certain  threshold  value  (2/3  of  the  amplitude  of 
each  variable).  Through  this  choice  of  simulated  activ¬ 
ity,  a  best  reflects  the  phase  difference  between  E  and 
M. 

Evening  and  morning  oscillators  are  identical  with  iden¬ 
tical  responses  to  light  pulses .  Differences  in  the  proper¬ 
ties  of  evening  and  morning  oscillators  in  hamsters 
have  been  suggested  by  the  differential  responses 
shown  by  the  beginning  and  end  of  the  activity  bands 
to  light  pulses  (Elliot  and  Tamarkin,  1994).  In  Figures 
7B  and  8B,  we  showed  that  these  responses  could  be 
simulated  well  by  two  identical  oscillators.  Due  to  the 
phase  difference  between  them,  light  hits  the  oscilla¬ 
tors  at  two  different  phases,  thus  evoking  distinct 
phase  shifts.  Although  a  more  complex  system  of 
oscillators  might  yield  similar  results,  we  have  shown 
here  that  a  simple  system  of  weakly  coupled,  identical 
oscillators  is  sufficient  to  explain  these  phenomena. 
We  further  assumed  that  the  tau  mutation  affected 
equally  the  periods  of  E  and  M  oscillators. 

CEMand  CMEhave  identical  absolute  values  and  opposite 
sign.  Asymmetry  either  in  the  oscillator  properties  or 
in  the  coupling  pathways  is  necessary  to  generate  a 
coupled  system  with  a  stable  phase  relationship  dif¬ 
ferent  from  0°  or  180°.  The  choice  of  two  identical  oscil¬ 
lators  coupled  by  CMEand  C£Mwith  the  same  absolute 
values  and  opposite  sign  (+,-)  was  again  chosen  for 
simplicity.  This  asymmetry  in  the  interactions,  one 
excitatory  and  the  other  inhibitory,  leads  to  some  par¬ 
ticular  features  in  the  model,  such  as  opposite  changes 
in  the  amplitudes  AE  and  AM  when  coupling  strength 
is  modified  (Fig.  5A,B).  If  the  signs  of  coupling  were 
equal— namely,  (+,+)  or  (-,-)— then  both  amplitudes 
would  change  in  the  same  direction  when  coupling 
strength  was  modified.  In  all  cases,  the  PRC  is  nega¬ 
tively  correlated  to  the  total  amplitude  of  the  system, 
which  is  the  sum  of  the  amplitudes  of  E  and  M  oscilla¬ 
tors.  This  is  shown  unequivocally  when  the  ampli¬ 
tudes  of  the  two  oscillators  change  in  opposite  direc¬ 
tions  (Fig.  5A,B). 

Aftereffects  of  a  light-dark  cycle  are  simulated  by  a  grad¬ 
ual  weakening  of  the  coupling  strength.  Aftereffects  are 
long-lasting,  slowly  decaying  alterations  in  pace- 
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maker  properties  caused  by  prior  light  schedules 
(Pittendrigh  and  Daan,  1976a).  The  phase  relationship 
between  E  and  M  depends  on  their  intrinsic  periods 
and  on  coupling  strength.  We  assume  that  the  strength 
of  coupling  is  high  in  hamsters  exposed  to  LD  cycles. 
When  the  condition  is  changed  to  DD,  the  coupling 
strength  slowly  decays  to  a  low  value.  Weakening  of 
coupling  leads  to  an  increase  in  the  phase  difference 
between  E  and  M,  reflected  in  the  increase  of  a.  We 
chose  not  to  assume  changes  in  oscillator  periods  dur¬ 
ing  prolonged  exposure  to  DD,  which  might  also  gen¬ 
erate  changes  in  a,  because  many  more  assumptions 
would  be  required.  For  the  purpose  of  this  analysis, 
our  simplest  assumption  of  coupling  strength  change 
alone  was  sufficient  to  explain  the  interconnection 
between  t,  a,  and  PRC  amplitude.  Our  primary  result 
that  system  amplitude  defines  the  PRC  amplitude  is 
independent  of  the  specific  model  assumptions  (Fig.^). 

The  Two  Views  of  Coupled  Oscillators 

Most  previous  modeling  studies  on  coupled  oscil¬ 
lator  systems  assumed  weak  coupling  and  strongly 
attracting  limit  cycles,  using  "phase"  oscillators  that 
lack  amplitude  information  (Daan  and  Berde,  1978; 
Kawato,  1985;  Mori  et  al.,  1994).  This  is  a  good  approx¬ 
imation  when  the  effects  of  coupling  on  the  amplitude 
of  the  oscillations  can  be  ignored.  Amplitude  can 
indeed  be  neglected  in  the  description  of  the  dynamics 
of  an  oscillator  if  it  remains  unchanged  by  experimen¬ 
tal  manipulations.  Amplitude  can  be  ignored,  in  par¬ 
ticular,  in  the  case  of  strongly  attracting  limit  cycles, 
which  are  able  to  restore  the  original  amplitudes 
quickly  when  perturbed  by  brief  external  pulses 
(Ding,  1987)  or  extremely  weakly  coupled  oscillators 
(Kuramoto,  1984). 

We  presented  two  different  views  of  coupled  oscil¬ 
lators,  the  system  view  and  the  distinct  oscillators 
view.  The  strength  of  coupling,  reflected  in  these  two 
views,  has  implications  for  (1)  the  degree  of  the  preser¬ 
vation  of  each  oscillator's  intrinsic  properties,  such  as 
free-running  period,  wave  shape,  amplitude,  and 
PRC,  and  (2)  the  degree  to  which  each  oscillator  moves 
independently  when  submitted  to  light  pulses. 

The  system  view  is  the  more  general  one,  in  which 
the  emergent  properties  of  the  system  as  a  whole  are 
characterized.  Amplitude  changes  due  to  strong  cou¬ 
pling  cannot  be  neglected  because  they  play  primary 
roles  in  the  behavior  of  the  system.  In  the  limit  of  weak 
coupling,  the  amplitudes  as  well  as  other  intrinsic 


properties  of  the  components  are  assumed  unaltered. 
In  this  limit,  we  can  see  immediate,  separate  responses 
of  distinct  oscillators  when  a  light  pulse  is  given. 
Moreover,  we  can  estimate  the  PRC  of  the  system  from 
the  intrinsic  PRCs  of  component  oscillators  and  repre¬ 
sent  these  components  by  phase  oscillators. 

The  pacemaker  in  hamsters  appears  to  lie  at  the 
border  between  weak  and  strong  coupling,  depend¬ 
ing  on  the  light  regime.  On  one  hand,  weak  coupling  is 
indicated  by  the  relatively  independent  phase  shifts  of 
the  component  oscillators  to  light  pulses  under  pro¬ 
longed  DD.  On  the  other  hand,  the  decrease  in  ATafter 
transferring  from  LD  and  during  long-term  exposure 
to  DD  (inferred  from  the  increase  in  PRC  amplitude) 
demonstrates  that  the  strength  of  the  coupling  cannot 
be  ignored.  No  sharp  line  divides  the  two  views. 
Either  view  is  useful  if  the  coupling  is  weak.  The 
choice  depends  on  the  specific  question  that  is 
addressed. 

Although  the  total  system  amplitude  does  not 
change  much  under  prolonged  exposure  to  DD,  espe¬ 
cially  for  the  wild-type  hamster,  the  amplitude  of  each 
component  oscillator  is  changed  (Fig.  5A).  Thus,  care 
has  to  be  taken  when  trying  to  measure  the  intrinsic 
properties  of  E  and  M,  such  as  PRC£,  PRCM,  xt.,  and  xM, 
during  transient  dissociation  of  activity  bands  and 
especially  when  splitting  occurs.  Although  each  oscil¬ 
lator  appears  to  be  independent  under  these  circum¬ 
stances,  each  is  in  fact  altered  by  the  effect  of  coupling. 
Therefore,  properties  measured  in  the  coupled  system 
may  not  reflect  precisely  those  intrinsic  to  the  oscilla¬ 
tors.  The  stronger  the  coupling  strength,  the  larger  are 
the  modifications  of  these  intrinsic  properties. 

Pittendrigh  et  al.  (1991)  proposed  that  the 
photoperiodic  time  measurement  in  Drosophila  would 
be  mediated  by  its  effect  on  pacemaker  amplitude. 
The  two-oscillator  system  provides  a  good  framework 
to  analyze  this  proposition  since  changes  in  the  phase 
relationship  between  E  and  M  clearly  occur  under  dif¬ 
ferent  photoperiods.  Elliot  (personal  communication 
and  abstract  from  the  1999  ICC  meeting)  has  per¬ 
formed  PRC  experiments  in  wild-type  hamsters 
exposed  to  different  skeleton  photoperiods.  Hamsters 
entrained  by  short  skeleton  photoperiods  have  dis¬ 
played  significant  increases  in  x,  a,  and  phase  shift 
responses.  These  findings  support  the  assumption 
that  different  light  regimes  can  affect  coupling 
strength  (which  define  the  phase  relationship  between 
£  and  M)  and  the  system  total  amplitude.  System 
amplitude,  in  turn,  affects  the  PRC  amplitude. 
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Period  Mutant  Pacemakers 

Because  x  and  PRC  define  the  phase  relationship 
between  the  clock  and  the  environmental  cycle  and 
because  this  phase  relationship  has  adaptive  value, 
Daan  and  Pittendrigh  (1976)  investigated  whether  a 
compensatory  mechanism  exists  between  the  two. 
They  found  indications  of  an  interdependence  of  land 
PRC  shape  that  provides  for  a  compensation  for 
day-to-day  instability  of  the  frequency  and  the 
interindividual  variation  of  x  within  the  rodent  spe¬ 
cies  they  analyzed.  This  compensatory  mechanism 
has  evolved  through  natural  selection  over  many  gen¬ 
erations.  Of  course,  a  period  mutant  organism  gener¬ 
ated  in  the  laboratory  is  not  adapted  to  a  specific  envi¬ 
ronmental  niche,  except  perhaps  by  chance.  On  the 
other  hand,  we  have  shown  that  physical  constraints 
exist  among  nonlinear  oscillator  properties,  and  these 
necessarily  have  to  be  obeyed  by  a  period  mutant.  In 
this  sense,  if  the  mutation  affects  such  biochemical 
parameters  as  synthesis  or  degradation  rates,  itneces- 
sarily  changes  both  the  period  and  the  pacemaker 
amplitude,  which  defines  the  PRC  amplitude.  Thus, 
we  have  pressure  exerted  by  natural  selection  for  a 
compensatory  mechanism  between  xand  PRC  as  well 
as  physical  constraints  defining  the  interconnection 
between  xand  PRC.  The  resulting  mechanism  reflects 
a  compromise  between  these  two  forces;  furthermore, 
this  compromise  can  be  modulated  by  a  separate 
mechanism  that  changes  the  sensitivity  of  the  system 
to  light  signals. 

Model  Predictions  and  Conclusions 

The  verisimilitude  of  computer  simulations  to 
physiological  results  can  never  prove  the  correctness 
of  a  particular  model.  Nevertheless,  the  close  agree¬ 
ment  between  our  simulations  and  the  available  data 
on  hamsters  leads  us  to  several  proposals.  First,  we 
strongly  support  the  proposition— namely,  that  the 
pacemaker  system  in  hamsters  comprises  two  oscilla¬ 
tors.  Second,  we  propose  that  the  amplitudes  of  E  and 
M  oscillators  in  tau  mutant  hamsters  are  smaller  than 
in  the  wild-type  animals.  Third,  we  propose  that 
weakening  of  coupling  between  E  and  M  during  con¬ 
tinuous  exposure  of  hamsters  to  DD  is  a  feasible 
hypothesis.  Finally,  we  argue  that  although  we  used 


equations  devised  for  the  Drosophila  pacemaker,  our 
analyses  do  not  depend  on  the  specific  properties  of 
any  system. 


APPENDIX 


We  consider  the  amplitude  of  each  oscillator  as  the 
sum  of  the  amplitudes  of  both  R  and  S  state  variables 
due  to  the  way  the  effect  of  light  pulse  is  incorporated 
into  our  equations.  In  some  models,  the  effect  of  light 
is  mimicked  by  an  instantaneous  change  in  the  level  of 
a  specific  state  variable.  In  that  case,  only  the  ampli¬ 
tude  of  that  variable  needs  to  be  considered.  In  our 
model,  we  represent  the  light  effect  by  a  1-h  change  in 
the  value  of  a  parameter,  affecting  all  state  variables, 
and  consequently  we  consider  all  their  amplitudes. 

The  simulation  change  between  wild-type  and  tau 
mutant  hamsters  was  made  by  altering  a  single-sys¬ 
tem  parameter.  Wild-type  oscillators  are  represented 
by  the  parameter  set  (a  -  0.85,  b  =  0.3,  c  =  0.8,  d  =  0.5) 
and  tau  mutant  oscillators  by  the  set  (a  =  0.85,  b  -  0.3,  c 
=  1.5 ,d-  0.5).  In  Figure  Al,  we  show  how  this  parame¬ 
ter  change  affects  x  and  the  amplitude  of  the  system. 
Each  line  in  the  figure  represents  the  new  values 
assumed  by  x  and  amplitude  when  one  parameter  is 
changed,  leaving  all  the  others  unaltered.  For  exam¬ 
ple,  the  point  indicated  by  "d  =  1.0"  represents  x  and 
amplitude  values  of  the  system  when  the  parameter 
set  is  (fl  =  0.85,  b  =  0.3,  c  -  0.8,  d  =  1.0).  The  wild-type 
hamster  is  represented  by  the  right-most  star.  Simi¬ 
larly,  x  and  amplitude  values  of  the  parameter  set 
defining  the  tau  mutant  system  are  indicated  by  the 
other  star.  A  change  in  parameter  b  in  the  range 
between  0.5  and  0.7  could  also  have  generated  a  sys¬ 
tem  with  x  =  20  h,  but  the  system  amplitude  would 
have  been  too  large  to  provide  a  type  0  PRC.  Two  main 
features  shown  by  this  figure  are  the  following:  (1) 
parameter  changes  that  modify  x  always  modify  the 
amplitude  of  this  nonlinear  oscillator,  and  (2) 
although  most  of  the  PRC  experiments  on  period 
mutant  organisms  have  shown  that  mutations  that 
shorten  x  also  decrease  amplitude  (since  they  show 
increasing  PRC  amplitude),  it  is  also  possible  to 
increase  amplitude,  as  shown  by  the  modifications 
generated  by  parameter  d. 
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Figure  Alinterdependenceoftand  amplitude  of  coupled  Pittend  righ-Pavlid  is  oscillators.  Each  parameter  set  (a,  b,  c,  and  d )  defines  a  spe¬ 
cific  land  amplitude  configuration  for  the  system.  The  wild-type  set  U  =  0.85,  b  =  0.3,  c  =  0.8,  d  =  0.5)  generates  a  system  with  T  =  24  h,  indi¬ 
cated  by  the  right-most  asterisk. The  tau  mutant  set  (a  =  0.85,  b  =  0.3,  c  =  1.5, d  =  0.5)  generates  a  new  configuration  with  t  =  20  h  and  smaller 
amplitude,  indicated  by  the  left-most  asterisk.  Each  parameter  is  changed  independently;  configurations  are  shown  with  respective 
parameter  values.  CMI  =  0.2  and  CEM  =  -0.2. 
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